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DEVELOPMENT OF AN ANALYSIS FOR 
THE DETERMINATION OF COUPLED HELICOPTER 
ROTOR/CONTROL SYSTEM DYNAMIC RESPONSE 

PART II - PROGRAM LISTING * 


By Lawrence R. Sutton 


Rochester Applied Science Associates, Inc. 

SUMMARY 


The computer program listing is presented for the program 
of investigating the effects of an anisotropically mounted flexi- 
ble swashplate on a deformed free-wake on the dynamic steady-state 
forced response characteristics of helicopter rotor systems. For 
the free-wake effects to be included, this program is used in 
conjunction with the free-wake geometry program listed in NASA 
CR-2111. The listing in this volume corresponds to the calcula- 
tions discussed in PART I. 


INTRODUCTION 


The determination of the effects of nonuniform control 
system (swashplate system or gyroscopic control system) along 
with the effects of nonuniform inflow due to free-wake on blade 
air loads and response for helicopters in steady maneuvers is 
essential in order to obtain more complete model of both aero- 
dynamics and elastomechanics of the helicopter rotor system. The 
program listed here includes the effects of an anisotropically 
mounted flexible swashplate (or nonuniform gyroscopic control 
system) and a deformed free-wake on the dynamic steady-state 
forced response characteristics of helicopter rotor systems. 

Three steps are necessary in obtaining blade response 
results including the effects of nonuniform inflow due to a 
free-wake and the effects of the nonuniformly supported control 
system. 

1. Preliminary calculations (or measured data) are used to 
define rotor system performance parameters and flight conditions 
Definitions of model parameters and program contol variables are 
necessary for program operation. 


*PART I - ANALYSIS AND APPLICATION is contained in NASA CR-2452. 



2. The results of the free-wake analysis are necessary to 
account for the effect of a nonuniform downwash distribution 
acting on the rotor system. For this the program listed in 
NASA CR-2111 is used to obtain wake-induced velocity influence 
coefficients (stored on BDSIG file) and bound circulations (stored 
on BDGAM file) for transfer to and use by this blade-swashplate 
dynamic response program. 

3. The blade response program listed in this report is then 
used to determine specified harmonics of blade flapwise, chord- 
wise, and pitching moments; blade deflections and twist, and 
their slopes; and blade spear and axial forces. This program 
also provides the corresponding harmonics of the swashplate 
response. Since the blade system responds to a downwash field, 
an iterative procedure is necessary to obtain the final dynamic 
response. 

Program input and output can be in English units or SI 
units. 
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OVERLAY (NS6 ,i.s,0) 

PROGRAM NS6 (IliPOT,OU TP UT , PUNCH ,CNE, TWO, THESE, EIGHT, NINE, TAPEl =ONB, 

1 TAPE2=T WO ,T APE3 =1 HRIsE, T APE4= PU NCH, T APB5=INPUT, TAPE 6 -OUTPUT, 

2 TAP'S8=EIGHT, TAPB9 = NINS) 

C MAINLINE PROGRAM FOE FORCED RESPONSE 

INTEGER CY18A 
INTEGER CY40,CY41>CY42 
INTEGER CY43,CY45,CY3,CYlS,CY44 
COMPLEX CYl , CY2 , CY 5, CY 6 
COMPLEX EPS (123) ,DBTSV 
REAL ITH(IO) 
c 

COM BON/AEF DB/DMS (4) 

COM MON/AER TP/AKT (4) ,ACT(4) ,AKP(4) ,ACP(4) 

COMMON/AKTAU/AKCI (4) ,TAU(4) ,SMLA{4) , AK(4) ,AC{4) ,CAPK ,CAPC 
COKHON/DIMS/NN1 ,MN2,NN3 
COMMON/BNAME/CS (4,20) ,SN (4, 20) 

COM MQN/RNAME 1/CS1 {4., 6) ,3Nl (4,6) 

COMMON/SCTA 3/CPSY (24,20) ,SPSY (24,20) 

COM MON /EP S A /BPS, DBYSV 

COMMON /CYI./CY3 ,CY!8, C Y 4 0 , C Y 4 1 , C Y 4 2 , C Y 4 3 , C Y 4 4 , C Y 4 5 , C Y 1 8 A 
COM MON /C YC/CY1 , C ?2 , C Y5, C Y6 

COM MC N/CYR/CY4 ,CY7,CY8, CY9 , CY19 , CYl 1 , CYl 2 ,C Yl 3, CY14 ,C Yl 5,CY1 6 , 

1 CY17,CY19,CY2&,CY21,CY?2,CY23,CY24,CY25,CY26,CY27,CY28,CY29, 

2 C Y30 , CY31 , C Y32, C Y33 , CY 34 , CY 35, CY36 , CY3 7 , C Y 38 , C Y.3 9, C Y46 ,C Y47 , C Y4 8 , 

3 CY4 Q ,CY51,CY5i,CY52,CY53 
CCMM0N/N01 / MXCPK , MXKQ , M XS11 E, MFASB, MXFA3 
CCM M0N/N02/MXCSB ,NFSBC,MXCFM, NO BC 
COMMON /N03/NS ,NSI ZEY , NFEA, NFS, MAXN ,NFPl 

COM M0N/N04/N COLS, N B, N F, S3IF :,N3 ISC, MXTKN, NIG 
CCMM0N/N05/NP, NSP, MODE, MFP, NOUT,MXQ,NAS ,N8S,NET,NPS 
COM MCN /NQ6/NEB D, NR IFC , MXSMI , NFL & F , MFLAP , NCT , MX rn 2Pl , NCS B 
COMMON /N 08 /JY HO 

COM MON/GYR/GYM , G YK ,GYC ,G K8 , GCB , GKP , GC P, G P P , GIX , GI Y , GIZ 

CCHM0N/CYM/CY54,CY55 

COM MON/IPHl /NN FA Z , NN B3 

COM MO N /I SMB/ MS 

C 0 M M 0 N / 1 B S I G / M I N P N , H A X P N 

CGMMCN/SW ASH/S WG J, SWSI ,SWM , SHE 

CO B MO N /V ELI / VK L (1C ,24) 

COM M0N/WA1 /CCN VG 
CCMMON/SA2./GAM (240) 

CCM MON/W A4/ALPH1 (10,24) , VEL F ( 1C , 24) , LTH 
COMMON/IMAI/IWAKB, NAFRO, MACRO 
CCM MON/ WA5/ ALL I 
C0MM0N/IWA3/MXIT 

. C 0 r MON/IW A2/ITI, N A S R , I MAX, X T £S T 
COM MON/IAR1/IAB (10) 

COM MO N /I W A4/NN4 , NN 5 , N N6 
COM M0N/SI1/SI (10) 
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CO ft MON/SXI3S/Y1 ,Y2, Y3 ,Y4 , Y5 
CQMMCN/IH/IH3 4 
COM MON /ASH DP/E J (4) 

C0MM1H/WA3/VZA1 (13,33) , V Y A 1 (13,33) ,VZPl ( 1 3 , 3 3) , V YP1 (1 3 , 33) 
COM M 0 N / D G a D / N H W I , N C V F , N P E S 
CCtitn N /VMAK2/V MAXI 
COMMON /IS 1/ N F F 
COt: MON/CYM .j /C.Y5.6 , CY57 
COP MO.S/HGAT/NWFC ,8P6 , STRT 
COM MO N/L I N A/I I N A , L A F I 
COM MON/C YM 3/C Y 58 
C 

100 IT 1=0 
NCUT =6 
MOQE=0 

1)0 1 11 = 1,10 

1 i as (id =0 

BEAD (5 , 941 ) H Hi ,NM 2 , S N 3 , N N 4 , N N 5 , N N .6 , N $ F C , N \i 'A. I , N C V F , H P S 5 
R E A P ( 5 , 901) Y1 ,Y2,Y3,Y4,Y5 

READ (5,941) I WAKE, NAE RO , H X IT , N .A S , KT ES T , I H 3 4 
IF ( IWAKE.I2Q . 0) GO TO 2 
WRITE (6,905) 

2 READ (5, 901) ALLI 

READ (5, 941) (IAP (II) , 11=1 , NA.ER0) 

»A».U = NAS*NAERO 
1 3 A X = 1 0 
DO 6 12=1,10 
6 SI (12) =0.0 

IF (I WAKE , »K,0) CALL OVERLAY (3LNS6 ,1 ,3 ,6H3LO&LL) 

CALL OVERLAY (3LNS6 ,1 ,1 , 6B RECALL) 

IF (MAKRO.EQ. t ) GO TO 10 
NORNH =NFF~ M F 

IF (MAXPN, LE. NOSSH) GO TO ID 

MAXPN=NCRNH 

WRITE (6,906) MAXPN 

10 IF (NWHI.LE. MAXPN') GO TO 11 
NHWI= MAXPN 

WRITE ( 6 , 907) NHWI 

11 CONTI HUE 

IF (IWAKE.EO.Ci) GO TO 66 
IF(NP.EQ.O) GO TO 12 
NP=0 

WRITE (6,918) NP 

12 CONTINUE 

IF(MCVF.GT.C) GO TO 20 
ncv F=ie 0 
W HI TS (6 , 9C 4 ) 

20 CONTINUE 

IF (ITI.EQ.O) GO TO 6 5 
74 CONTINUE 
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IF ( K A m 0 , i; 0.0) G 0 TO 66 
DO 7 17-1 ,1*. 

7 S I ( 17) =C . f 
DO 3 MS=1,NB 

3 CALL AEBQB 
BEHIND 2 
REWIND 3 

66 MFP=~ 

IF(NP.GT, MAXPN) GO TO 5)0 
IF (I W A KB. I! Q. ij ) WHITE (6,940) UP 
CYl =CY2*NP*C Y4 
IF(NB.GT.I) GO TO £5 
MS-1 

CALL OVERLAY (3LNS6,1 , 4 , 6HRBCALL) 
CALL OVERLAY (3LNS6 ,1 , 5, 6HREC ALL) 
CALL OVERLAY (3LNS6,1 , 6 , 6HEFC ALL) 
CALL OVERLAY (3LNS6 ,2,0 , 6HRBCALL) 
MFP=1 

IF (IHAKE.EQ.G) GO TO 4 
IF (N PBS. HO. 0 ) GO TO 13 
WHITE (6, 94(f) NP 
GC TO 4 

13 I F ( IT I . GK . M X IT ) WRITE (6 ,940) NP 

4 CONTINUE 
REWIND 2 
REWIND 3 

CALL OVERLAY ( 3LNS6 ,1 , 4 , 6HRECALL) 
GO TO 69 

65 DO 67 MS=1,NB 

CALL OVERLAY {3L.N56 ,1 ,4, 6HSIC&LI-) 

67 CALL OVERLAY (3LNS6 ,1 , 5, 6HRKCALL) 
CALL OVERLAY. (3LNS6 ,1,6, 6HSECALL) 
CALL OVERLAY (3LNS 6, 2,0, 6HR EC A LL) 
MFP=1 

IF (IWAKE.FQ.fi) GO TO 5 
IF(NPRS.EO.O) GO TO 14 
UNITE (6# 940) NP 
GC TO 5 

14 IF (ITI. GF.MXIT) WE ITH (6 , 940) NP 

5 CONTINUE 
REWIND 2 
REWIND 3 

DO 68 MS=1 , NB 

68 CALL OVERLAY (3LNS6 ,1 ,4, 6HBF CALL) 

69 IF ( IHAKE.BQ, 0) GO TO 82 
IFfNP.LT.NHWI) GO TO 02 
IF (ITI.EQ.MXIT) GO TO 82 

CALL OVERLAY (3LNS6/1 , 2 , 6HREC ALL) 

NP=G 

GO TO 81 
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82 NF=NP+1 
REWIND 2 
REWIND 3 
GO TO 66 
81 CONTINUE 
REWIND 2 
REWIND 3 

IF (CONVG.Gl.ALLI) GO TO 95 
ITI=MXIT 
GO TO 98 
95 I TI =1 TI + 1 

98 WRITE (6,993) ITI, ALLI ,CCNVG 

CALL OVERLAY (3LNS6 ,1 , 3, 6HRICALL) 

GO TO 7a 

500 READ (5,941) IGO 

IF (IGO.EQ.f) STOP 
GO TO 100 

901 FORMAT ( 8 F10.0) 

905 FORMAT (/36X, *IF IWAKE EQUALS 1, SIGMA AND GAMMA AFP? Y 5 MUST BE ATT 
1 ACHED*) 

904 FORMAT {/36X, ^SOLUTION WILL NOT CONVERGE WITH NCVF = 0, NCVF RESET 
1 = 100 *) 

906 FORMAT (/25X , AXPN HAS BEEN RESET TO THE MAXIMUM VALUE ALLOWED BY 
1 THE VALUE OF N FF , MAX E N=* , 13) 

507 FORMAT (/4lX,*NHWI CANNOT BE GEE ATI. E THAN MAXPN.NHWI RESBT=*,I3) 

908 FORMAT (/40 X , *FOR WAKE ITERATION NP MUST START AT C,NP FES?T=*,I3) 

941 FORMAT (1615) 

940 FORMAT ( 1 Hi , 56X, 2 OH FORCED RESPONSE, NP=, 12 ,/59X, 13HSTATF VECTORS/) 

993 FORMAT (/33X , 5HITI =, 13 , 10 X , 6HALLI =,E15.7, 5X, 7HERROF =,E15.7) 

END 
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SUBROUTINE AEEOM 
REAL TC (34) 

REAL TMO (24) 

COMPLEX AMF (1230) 

COMPLEX AHA (1700) 

COMPLEX CY19A 

COMPLEX CY1 ,CY2,CY5,CY6 

INTEGER CYl 8 A 

INTEGER CY 40,€Y41,CY42 

INTEGER CY43,CY45,CY3,CY18, CY44 

REAL LTH(IO) 

C COMMON AREAS 

C 

COM MON/C YI/CY3 , C Yi 8 , €¥ 40 , CY4 1 , CY42 , CY 4 3 , CY4 4 , CY 45 ,C Yl 8A 
CCM MON/C YC/CY1 ,CY2,CY5,CY6 

COM MON/C YR/CY4 , CY7 , CY8 , CY9 , CY 10 , CYl 1 , CYl 2, CYl 3 , C Y 1 4 , CYl 5 , CYl 6 , 

1 CY17,CY19,CY20,CY21,CY22,CY23,CY24,CY25,CY26,CY27,CY28,CY29, 

2 CY3 f 3 , CY31 , C Y32, CY3 3 , CY34 , CY 35, CY36 , CY37 , CY38 , C Y 3 9, CY4 6 ,CY47,CY48, 

3 CY 4 9 , C Y 5i) , C Y 5 1 , C Y 52 , C Y 5 3 
COMMON/NOl /MXCPK,MXKQ,MXSMB, MFASB, MXFAB 
C0MM0N/NQ2/MXCSB ,NES8C , MXCPM , NEBC 
COMMON /NO 3 /NS ,NSIZEY ,NFE'A, NES , HAXJf, NFP 1 
COMM 0N/N04/NC0LS ,NB ,NF ,NEIFC , HE ISC , MXTKN , NIG 

C o M H 0 N / N 0 5 / N P , N S P , H 0 D E , M F P , N 0 U T , M X Q , N A S , N 3 S , N 1 T , N P S 
COM M 0 N./N 06 /N R B D , NR IF C , M X S M I , NF L A P , MF L A P , NCT , M XT 2 P 1 , N CS B 
CCMM0N/I3RIG/MINPN, MAXPN 
C0MM0N/WA4/ALP Hi (10,24) , VELF (TO ,24) , LTH 
CCMMCN/VBL 1 /VEL (13,24) 

COMMON /IS1/NFF 

COM MON/DI MS /N N 1 , UN 2 , N N 3 

COM MON/SCTAB/CPSY (24 ,20) ,SP3Y (24,20) 

COMMON/SUB/Y (225) 

COM MON /S A IN/SP (742 5) 

COM MON/A MFT/AMF 

COMMON /AM AT/AM A 

CCMM0N/IWA1/IWAKE, NAERO, MACRO 

COMMON/C YM /€ Y54,CY55 

COM MO N/IWA4/N N4 , NN 5, N N6 

COM MON /ISM 3/MS 

CCMM0N/IWA7/ITI, NAN R , IM AX, KTSS7 

COMM 0 8 /I L 1 /I L 

COM MON/C YMl/CY 56 , CY57 

CO M MO N /N G A T/N W FC , N P 6 , N T R T 

COM MON/C YM 2/C Y 58 

IF(IMAKE.EQ.O) GO TO 4 

I F ( MS . GT, 1 ) GO TO 4 

DO 3 15 = 1 , NAERO 

DO 3 J5=1 , NAS • 

3 VEL (15, J5) =CY4*CY50*VEL (15, J5) 

4 CONTINUE 
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DO B v75=1 , NAS 

8 TMO (J5) =0,0 
MXNFF=2*NFF+1 
IL=C» 

IF (NET . F 0 . C ) GO TO 9 

BEAD (5,703) ( (VEL (IRITB , JPIT2) , IEIT 1=1 , SAERO) , JHIT1=1 ,NAS) 

9 IF ( NN6 .BQ, 0) GO TO 810 

WRITE (6, 801 ) ( (VBL (1 1 ,Kl) ,I1=1,NAEB0) ,Kl=1,NAS) 

801 FORMAT (* INDUCED VELOCITIES FROM ABFCM */ (5016.7) ) 

810 DO 1j 1=1,34 
10 TC (I) = 0. 

I F (I W A KB . EQ. 0 } GO TO 2 
I F ( MS . GT . 1 ) GO TO 6 
2 DO 12 I = 1 , N N 1! 

12 AM A (I) = CMPLX (0. ,0. ) 

6 DO 5 1=1, NN2 

5 AM? (I) =CMPLX (C*. ,0.) 

IF (ITI.EQ.G) GO TO 7 
HEAD (2) SD 

7 DO 26 1=1, NS 
L=(I-1) * NSIZEY 

Y (203 ) = SD (L+20,3) 

IF (Y (203) -1) 26,13,26 

13 I L=IL+1 

14 DO 16 J = 1 , N SI7.FY 

16 Y (J) =SD (L+J) 

LTH (IT ) = Y (7 5) 

DC 24 K=1 , NAS 

IF (I WAKE . EQ. 0) GO TO 17 

Y (1 80 ) = VEL (II , K) 

17 CONTINUE 

IF (NET . EQ. 0) GO TO 18 

Y (180) =VBL{ IL,K) 

18 CONTINUE 

CYl 9= (K-1) * (2. 0*.CY12) /NAS 
CY20 = SPSY (K, 1) 

CY21 = CPSY (K, 1) 

CY53=-CY47* (CY21 Y51 +C Y20'*CY52) 

SC Y53 =SI N (C Y53 ) 

C CY 5 3 =CO S ( C Y 5 3 ) 

SANG = Y (19) -*CCY53+Y (20) *SCY53 

CANG=Y (20) ^CCY53- Y (19) *SCY53 

PTJLL=CY54+CY55*CY56 

R0LL=CY4 8*CY1 6-CY5 5’* t CY57*CY 1 7 

TrjRN=CY48*CYl7+CY55*€Y57*CYl6+CY4 

CY22 = CYl 4 * CYl 6 * C Y >J 

CY23 = CYl 4 * C Y 1 6 * CY21 

CY2 4 = CYl 4 * CYl 7 

CY3 5 = Y (23) * SAN G- Y (11) *C ANG 

C Y36 = Y (^2) *CANG+Y (24) *S ANG 
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CY37=y (11) *SAKG + Y (28) *CAHG 

CY38=Y (24) *CANG-Y (12) *SAKG 

CCS =P!JLL*CY21 + S0LL*CY20 

CCS=RQLL*CY2 1-PULL *C Y20 

C Y 2 5 - T U F< N * Y ( 1 6 ) * S A N G + C C A * C Y 3 5 + C C B *' C Y 3 6 

CY26=TURN*Y (1 6) *CA NG + CCa tt C Y37+CC3*C Y38 

C Y27=CY? 2 

CY28=CY23+TUBN*CY58 
CY29 - Y (78) * CY25 
CY30 = Y (78) * C Y ? 6 
CY31=CY24+Y (1 80) -CC3*CY58 
CY32=CY31*Y (16) #SA NG 
CY33=CY3 1*Y (16)*CANG 
CY34=CY31*Y (15) 

CY39=TURH*Y (1 5) -CCA*Y (29) -CCB*Y (25) 
SSS1=CY27*CY35+CY28*CY36+CY32+CY?6* (Y { 3) -3Y58) 
SSS2=CY27*CY37+CY28*CY38+CY33-CY39*Y (4) -CY25* (Y (3) -CY58) 
SSS3=CY27*Y (29) +CY28*Y (25) -CY34»CY26*Y (4) 

S3S4=S3S3+CY3C 

Y (8 2) --CY39 

Y (83) =3551 

Y (84) =SSS2+ Y (78) *CY39 
CALL ASK 0 

IF ( NN6 . BQ. 0 ) GO TO 805 

WRITE (6 , 802) Y (83) ,Y{84) ,Y (90) , Y{180) ,CY14 

802 FORMAT {* BLADE VELOCITIES F30M AEROM */(5Gl6.7) ) 

W RITE (6 , 80 3 ) Y (1 24 ) , Y (1 2 5) , Y { 1 2 6) 

803 FORMAT (* BLADE STEADY LOADS PER UNIT LENGTH */(5Gl6.7)) 
805 IF (I WAKE. BO. 0) GO TO 113 

ALP HI (IL , K) =Y (91) 

VELF (IL , K) = Y (90) 

117 FORMAT ( 10Gl 2. 4) 

118 CONTINUE 

C CREATE MATRIX C, 

TC (1 ) = Y (123 ) * Y (78) - Y(117) 

TC (2) = -Y (120) 

TC (3) = Y (123) 

TC (4) = - Y (121 ) * Y (78) + Yfl.15) 

TC ( 5) = Y (1 18) 

TC (6) = -Y (121) 

TC (7) = Y (12 2) * Y (78) -Y(116) 

TC (8) = — Y ft 19) 

TC (9) = Y (122) 

TC (10) = Y (124) 

TC (11) = -Y (12 5) 

TC(12) =-Y (1 20) *C 12 6 + Y (123) *CY25 
TC (13) =Y (123) *SSS1-Y (120) *SSS2 
TC (14) = Y (123) * CY39 

TC ( 1 5) =+ Y (120) *SSS3 + Y(123) *C Y29-Y (117)*CY25 
TC (16 ) = Y (120) * CY 39 
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TC (17)=-Y (123) *SSS4 + Y (117) *CY26 
TC (18) = Y (126) 

TC (19) =+Y (118) *CY26-Y (121) *CY25 
TC (20) = Y (118)*SSS2-Y (121) *SSSl-Y (125) 

TC (21) = -Y (121) * CY39 

TC (22) =-Y (118) *SSS3-Y (121) *CY29+Y (115)*CY25 
TC (23 ) = -Y (118) * CY39 
TC (24) = + Y (121) *SSS4-Y(115) *CY26 
TC ( 2 5 ) = - Y (126) 

TC ( 26) =- Y (119) *CY26 + Y (122) *CY25 
TC (27) =Y (122) *SSS1-Y (119) *SSS2-Y (124) 

TC (28) = Y (122) * C Y 39 

TC (29) =+ Y (119) *SSS3+Y (122) *CY29-Y (116) *CY25 
TC (30) = Y (119) * CY39 

TC (31) =-Y(122)"*SSS4 + Y(116) *CY26 
TC (32) = -Y (126) 

TC (33 ) = Y (124) 

TC (34) = -Y (125) 

DC 22 NN-1 , MXSMI 
N=NN-N.FPl 

CALL PS YK (CPSY,SPSY,K,N,CY19A) 

IPO=34* (TL-1) *HXSMI+ (N'N-1) *34 
IF (IWAKE .EQ. '}) GO TO 51 
IF (MS . GT. 1) GO TO 22 

51 DC 20 J = 1 , 34 

A MA (I PQ+ J) = A Pi A (IPQ + J) - Y(7 5) * TC (J) * CY19A 
20 CONTINUE 
22 CONTINUE 

DO 220 NN=1,MXNFF 
N=NN-NFF-1 

CALL PSYK (CPSY,SPSY,K,N,CYl 9k) 

JPQ=3* (IL-1) *MXNFF+ (NN-1 ) *3 
DO 220 J =1 , 3 

220 AMF (JPQ+ J) = AMF (JPQ + J) -Y (75) *TC ( J + 31) +CY19A 
IF (NTBT, NB. 1 ) GO TO 24 

TMO (K) =TfiO(K) -Y (3) *Y (75) * (Y (124) *Y (35) +Y (1 2 5) *Y (36) ) 

24 CONTINUE 

26 CONTINUE 

IF ( NTET. NE. 1 ) GO TO 808 
WRITE (6, 807) (TMO (K) ,K = 1 ,NAS) 

807 FOB MAT {* BLADE THRUST MOMENT AT EACH AZIMUTH */(5Gl4.5) ) 

808 CONTINUE 
J=MXSMI*34*IL 

IF(IWAKE .EQ. 0) GO TO 52 
IF (MS .GT. 1) GO TO 53 

52 DO 28 1=1, J 

28 AMA(I) = AKA (I) / NAS 

53 JF=MXNFF*'3*IL 
DO 228 1=1, JF 

22 8 AMF (I ) = AMF (I) /NAS 
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IF (I WAKE .EQ. 0) GO TO 35 
WRITE (3) AMF 
GO TO 36 

35 WRITE (3) AMA , AMF 

36 CONTINUE 

9876 FORMAT (*0* , 10Gl 2 . 4) 

703 FOR MAT (8 FlO. 6) 

902 FORMAT (201,215, 3115, 7) 

903 FORMAT (10E1 3. 6) 

39 RETURN 

END 



SUEEOOTINE AERO 
CCMFLKX CYt ,CY2,CY5,CY6 
INTEGER CY43,CY45,CY3,CY18,CY44 
INTEGER CY40 , CY41 , CY42 
INTEGER CY18A 

COMMGN/CYI/CY3 , C Yl 8 , CY40 , CY41 , CY42 , CY 43 , CY4 4 ,CY45 ,C Yl 8 A 
COMMON /C YC /C Yl ,CY2,CY5,CY6 

COMMON/C YE/CY4 , CY7 , CY8 , CY9 , CY19 , CYl 1 , CY 1.2 ,€Yl 3 , C Y14 , CYt 5, CYl 6 , 

1 CY17,CY19,CY2€?,CY21 , C Y2 2 , C Y23 , CY24 ,C Y2 5 ,0726 , C Y27 , CY28 , C Y29 , 

2 CT30,CY31 ,CY32,CY33,CY34,CY35, CY36 , CY3 7 ,CY38 , C Y39, CY46 ,CY47,CY48 , 

3 C Y 4 9 , C Y 50 , C Y5 1 , C Y 52 , C Y 5 3 
COMMON /SUB/Y (225) 

C08H0H/IH/IH34 
COMMON/I II /IL 

COMMON /Sit /SI (10) 

COMMON/S UBS/Yl ,Y2,Y3,Y4,Y5 
EMT=CY4*CY50/CY9 

Y (90) = SQRT {Y (83) * Y (83) + Y(84) * Y(84)) 

UN= Y (90) / (CY5G*CY4) 

IF (Y (84) ) 222,112,111 

111 IF (Y (83) ) 555,556,666 
222 IF (Y (83) ) 333,445,444 

112 IF (Y (83) ) 113,114,114 

113 Y (91) = CY12 
GO TO 668 

114 Y (91) =0,0 
GO TO 668 

333 Y (9 1) = C Y 12“ AT AN { A B S ( Y (84) / Y (83) ) ) 

GO TO 668 

444 Y (91) = ATAN (ABS (Y (84) /Y (83) ) ) 

GO TO 668 

445 Y (9 1) = .5 *CY12 
GO TO 668 

555 Y (91) = CYl 2 + ATAN (ABS (Y(84) /Y (83) ) ) 

GO TO 668 

556 Y (91) =1. 5*CY12 
GO TO 668 

666 Y (91) = 2. *C Yl 2 - ATAN (ABS (Y (84) /Y (83) ) ) 

668 CONTINUE 

Y (92) = Y (90) / CY 9 

Y (93) = 1. / (1. - Y (92) * Y (92) ) 

Y (94) = (2, - Y (92) * Y(92)) * Y(93) 

CALL NS1214 (EMT , UN, Y (9 1 ) ,Y{95) ,Y(96) ,Y (99) , Y (97) ,Y(100) ,Y{98) ) 

IF (IH34.E0.0) GO TO 760 
IF (CY19.LT. 3.0) GO TO 760 
IF (CYl 9 . GT .4.8) GO TO 760 
AS= (. 22688 + Y5) * (1.0-Y (9 2) ) 

IF(Y (91) .LT.AS) GO TO 760 
IF (SI (IL) .NE.C.O) GO TO 20 
SI (IL) = ASIN ( (Y (95) -Y3) /Yl) /Y2 
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20 Y (S 5) =Y3 + Y1 *SIN (Y2* (CY19-SI (IL) ) ) 

Y (99) =.5'* (1.0-COS ( Y 2 * (CY-19-SI (II) ) ) )*Y4 

760 Y (101) = CI10 * Y { 7 9 ) / 2. 

Y (102) = Y (101 ) * Y (90) 

Y (103) = Y (101) / Y ( 90) 

Y (104) = Y (1 01) * Y (79) 

Y (105) = Y(104) * C Y 1 2 / 2. 

Y (1 06) = Y (79) / 4. - Y (78) 

Y (1 07} = Y (84) * Y (9 5) 

Y (108) = Y (84) * Y (97) 

Y (109) = Y (83) * Y (95) 

Y (1 10) = Y (83) * Y (97) 

Y (111) = Y (84) * Y (96) 

Y (1 12) = Y (84) * Y (98) 

Y (113) = Y ( 83) * Y (96) 

Y (114) = Y (83) * Y (98) 

Y (1 15) = -Y (10 5) * Y (84) 

Y (116) = Y (115) * ¥(8 3) 

Y (1 17) = -Y (78) * Y (116) 

Y (118) = -Y (1C 3) * ( Y (84) * ( Y (109) * ¥(93) + Y(111)+ Y (108) ) 

1 + Y ( 83) * ( Y (9 4) * Y (IIP) + Y (112) ) ) 

Y (119) = Y (103) * ( Y (83) * (Y{94) * Y (109) + Y (1 11 ) ) -Y (84) 

1 * ( Y (1 10) * Y (93 ) + Y (112) - Y (107) ) ) + Y (105) * Y(82) 

Y (1 20) = Y (104) * (Y (94) * Y (83) * Y(99) + Y(84) * Y (190) ) 

1+Y (106) *Y (119) -Y (105) *Y (82) (Y (106) +Y (78) ) 

Y (121) = -Y (103) * ( Y (84) * (Y (94) * Y (107) - Y (1 13) ) 

1 + Y ( 83) * { Y (1)8) * Y { 93) + Y (109) - Y(114))) -Y (105) * Y (8 2) 

Y (122) = Y ( 1 C 3 ) * ( Y (83) * ( Y (1G7) * Y(93) - Y(113) - Y (111) ) 

1 - Y(84) * ( Y (94) * Y (103) - Y <1 14))) 

Y (1 23) =Y (10 4) *■ (Y (94) *Y (84) *Y (9 9) -Y (83) *Y (190) ) 

1 + Y (106) * Y (122) 

Y (124) = -Y (102) * ( Y (107) + Y(110)> + Y(115) * Y(82) 

Y (125) = Y (102) * (Y (109) -Y (1C 8) ) + Y(116) * Y(82) 

Y (126) = Y (102) » Y ( 7 9) * Y (90) * Y (99) + Y(1)6) * Y (125) 

1 + Y (117) * Y (82) - Y (106) * Y{116) *Y (82) 

SCY53=SIN (CY53) 

CCY 53 =COS (CY53) 

AAl 1 5 = Y (115) 

AAl 1 6=Y (116) 

AAl 18-Y (118) 

AAl 1 9 =Y (119) 

A A1 21 = Y (121 ) 

AAl 22 = Y (122) 

A Al 24 = Y (124) 

AAl 25=Y (125) 

Y (1 1 5 ) =AAl 1 5*CCY53-A A1 1 6*SCY53 

Y (116) =AA116*CCY53+AA11 5*SCY53 

Y (118) =AAl18*CC'Y53-AAlt9*SCY53 

Y (119) =AA11 9*CCY53+AA118*SCY53 

Y (121) =AA121*CCY53-AA122*SCY53 
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Y (122) »AAl2.2*CCY53+AAl21*SCY53 

Y (124) «AA124*CCY53-AA125*SCY53 

Y (1 25) =AA125*CCY53+AA124*SCY53 
BE TURN 

670 FORMAT (1 HO ,* Y (91 ) = * ,E1 2. 4 ,3X,*Y (92) = *,El 2. 4 , 3X,*Y (95) = *,£12.4) 

671 FOR MAT (1 H ,* Y(97) = *,10E12.4) 

672 FORMAT (1 H ,* Y (9 9 ) = *,10JS12.4) 

END 
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SUBROUTINE NS1214 (EMT , U , APHIJ, CLIFT , ASLO P,€MOME , CDR AG, CMSL, CDSL) 
COMMQN/LINA/IINA ,LAFI 
C 

C SUBROUTINE TO COMPUTE CLIFT-LIFT COEFFICIENT 

C AS LOP = LI FT CURVE SLOPS 

C CMOME=MOMENT COEFFICIENT 

C CBRAG* DRAG COEFFICIENT 

C CMSL= MOMENT CURVE SLOPE 

C CBSL=DRAG CURVE SLOPS 

C FORMULAS TAKEN FROM CURVE FITS BY P.C. 

C 

CLIFT=0. 

A SLOP-0 • 

CMOME=0. 

CDR AG=0 • 

CMS L=0. 

CDS L~0 • 

C 

130 NEG =1 

EMIJ=EMT*ABS (U) 

SQT=SQ.RT (1 • — EMI J*E MI J) 

Cl =1. -EMIJ 
C2=. 22689*C1 
IF (IINA.FQ.O) GO TO 97 
C2*C2* (1.0 + . 01 *LAFI) 

97 IF (APHIJ) 181,18 2, 182 

181 APHIJ*- A PHI J 
HEG=-1*NEG 

182 IF{APHIJ-3. 1415926) 184,184,183 

183 APHXJ=APHIJ-3. 1415926*2. 

GO TO 97 

134 I F (APHI J-C2 ) 185,187,167 
185 ASLOP=5. 7296/SOT 
Cl I FT =A SLOP* APHI J 
CDR AG=. 006 + • 13131* APHIJ* APHIJ 
C MO MS *0.0 
CM$L=0..0 

CDS L=» 2 6262* APHI J 
GO TO 250 

187 IF (APHI J-. 34906) 189,191,191 
189 C L I F T = . 2 9 2 6 9 * C 1 + ( 1 . 3 * E M I J- , 5 9 ) * A P H I J 
C MO Mg =0.0 

C2= (. 12 2 17+. 2268 9*. EMIJ) *SQT 
CLI FT=CLI FT/C2 
A SLOP* (1 »3*EMIJ-. 59) /C2 
C MS L=C • 0 
GO TO 210 

191 IF (APHI J-2. 740'2) 193,195,195 
193 S-SIN (APHIJ) 

S2=SIN (2.*APHIJ) 
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S3-SIK ( 3,*APHIJ) 

S4 = SI N (4. *APHI J) 

CL I FT = { . 0 80 3 7 3* S + 1 . 0 4 30 8* S2- . 01 1 9 59 *S 3 + . 0 2 3 1 2 7* S 4 ) /SOT 
C=COS (APHI J) 

C2=CQS <2.* A PHI J) 

C3=COS(3.*APHIJ) 

€4= COS (4.* APHI J) 

ASLOP= { . 080373*0 +2.0861 5*C2- . 03 31 77*03+ . 09250 8*C4) /SOT 
CDBAG = (1.1233-.C2989 4*C-1.r 0 603*02 + .1)031 15* 03- . 09 148 7*C 4) /SQT 
CBSL = (• 0 298 9 4*S + 2..-01 206*52-. fC9345*S3+. 365948*34) /SQT 
CMOME = <-.02827*S+.14 022*S2-.30622*S3 + . 01012*34) /SQT 
1 -. 25* (CLIFT*C+CDPAG*S) 

CMSL= {-* 56986$*C+. 2733 5*C2- . 028178*03+ , 94246 6*04 + S* (. DC 9263* S 2 
1 - . 0 O 3 97 2* S3— * 0 8 54 5* S4 ) ) /SQT 
GO TO .250 

195 IF (APHIJ-3* 0020) 1 97,1 99,199 
197 C LI F T =- (.4704 + . 1031 3*APHIJ) /SQT 
ASLOP=— • 1031 3/SQT 
CMOME--. 3452984/SQT 
CMSL=0. 

GO TO 210 

199 IF ( APHIJ-3. 141 5926) 200,2*0,260 

200 CLIFT= (-17. 550 + 5. 5964*APflIJ)'/SQT 
ASLO? = 5. 5864/SQT 

CM0M3*-. 3452934* (1.0- ( (APHIJ-3. 002) /. 13 95927) ) /SQT 
CHSL=. 34529 84/ (.1395927* SQT) 

210 CDRAG = (1.1233-.C29894*COS(A?HIJ) -1. 306C 3* COS (2.*APHIJ) 

1 +.0031 1 5*COS (3.*APHIJ) -.C91487*COS (4.*APHIJ) ) /SQT 

CD SL =(,*29894*3 IN (APHI J) +2. 01206*SI N (2. *APHIJ) 

1 -.OC 9345*SIN (3 . * APHI J) +. 36 5943* SIN (4.* APHI J) ) /SQT 

250 IF(NBG) 25 5,255,266 
255 CLIFT-- CLIFT 
C MO ME =- C M.Q M E 
APHU =- APHI J 
CDSL=-CDSL 
260 CONTINUE 
C 

300 CONTINUE 
BSTUEN 
END 
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SUBROUTINE PSY K (CPSY , SPS Y, K , Nl , CYl 9) 
COMPLEX CYl 9 

DIMENSION CPSY (24,20) ,SPSY (24,20) 

N = IABS (Nl) 

IF (Nl.BQ.fl) GO. TO 110 
A = CPSY(K,N) 

B = SPSY(K,N) 

110 IF ( Nl ) 120,121,122 

120 CYl 9 = CBPLX (A ,B) 

GO TO 123 

121 CYl 9 = CMPLX {1 .0,0,0) 

GO TO 123 

122 CYl 9 = CMP LX (A,-B) 

123 CONTINUE 
EETURN 
END 



OVERLAY (NS6,1,G) 

PROGRAM DUMMY 
COMPLEX FKN (123) 

COMPLEX FAS B (25) ,CTSB (25) ,FLSB (25) ,FLSC (25) »FASD(25) ,CTSD (25) , 
1 FASC (25) 

COMPLEX B (1800.) ,SHLB (303) ,SMLC(300) , SMLD (300) 

COMPLEX FAB (150) , FLB (1 50) ,CTB (1 50) 

COMPLEX H (60) , CTH (5) , FLH (5) , FAH (5) 

COM MON/FKNl/FKN 

COM MON/BTSl /B r SMLB ,SMLC,SMLD 

COMMON /BTS 2/ F A B r FL B , CTB 

COMMON/BTS3/FASB,CTSB,FLSB, FLSC, FASD, CT S D f FASC 
COM MOH/RETS/H , CTH , FLH , FAH 
COM M0N/R07 /I , •! 

RETURN 

END 
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a n n 


OVERLAY (856,1,1) 

PROGRAM SETUP 

INTEGER CY18A 

INTEGER CY40,CY41,CY42 

INTEGER CY43,CY45,CY3,€Yl8,CY44 

COMPLEX CY1 ,CY2,CY5,CY6 

COMMON AREAS 


COM MON/SWASH/S WGJ, SWEI ,SWM ,SWR 
COM MO N/B NAM E/CS (4,20) , S N (4 , 2 0) 

COM MON/BN AMEl/CSl (4,6) ,SN1 (4,6) 

COM MON/SCTAB/CPSY (24,20) ,SPSY (24,20) 

CCMMON/SAIN/SD (7425) 

COMMON/SUB/Y (225) 

COM MON/AKTA U/AKCI (4) ,TAU(4) ,SMLA(4) ,AK{4) ,AC(4) ,CAPK,CAPC 
CCMMON/GYR/GYM ,GYK ,GYC ,GKB,GCB, GKP , GC P, GRP , GIX , GI Y , GIZ 
COMMON /C Yl / C Y 3 ,CYl 8 , CY40 ,CY4 1 , CY42 , CY43 , CY44 ,CY 45 ,C Yl 8A 
COM MON/C YC/CYl ,CY2,CY5,CY6 

COM MON/C YR/CY4 ,CY7 ,CY8 , CY9 , CY10 , CYl 1 , CY 1 2 , CY1 3 , CY1 4 ,C Yl 5,CY1 fe , 

1 C Yl 7 , CY 1 9 , C Y2f , C Y2 1 , C Y2 2 ,C Y23 , C Y24 , CY2 5 , CY2 6 , C Y27 , C Y28 , C Y29 , 

2 CY30,CY31 ,CY32,CY33,CY34,CY35,CY36,CY37,CY38,CY39,CY46,CY47 ,CY48, 

3 CY 49 , CY 50 ,CY51,CY52,CY53 
COM MON /C YM/CY5 4 , CY55 
COMMON/IWA1/IWAKE, NAERO, MAERO 
CCMMON/I WA2/ITI , NANR, IMAX, K.TESf 
CO M MO N /I W A 3/M X I T 

C C M M 0 N / 1 W A 4 / N N 4 , N N 5 , N N 6 
CO M MO N/DI MS /NH 1 , NN 2 , N N 3 
COM MO M/I PHI / NNFAZ , NNBS 
COM MON/I Si /NFF 

COMMQN/NOl /MXCPK ,MXK3 , MXSM E, MFASB,MXFAB 

COM MON/N02/MXCSB ,NESBC, MXCP M , NEBC 

COMMON /N03/HS ,NSI ZEY , NFi'A, MSS , MAXN , NFP 1 

COMM 0N/N04./N COLS , NB , N F , NEIFC , NEISC , MXTKN , NIG 

CCM M0N/N05/N P , NS P, MODE, MFP, NOUT , M.XQ , N AS , MBS, NET , NPS 

COM M0N/N06/NEBD, NF IFC , MXSM I , NFL A I, MFLAP, NCT , H XT 2 Pi , NCS B 

COM M0N/NO8/JYRO 

COM MO N/I BF; IG/M IN PN , MA XP N 

COM MO N /AER D M/D MS (4) 

COM MON/AERTP/AKT (4) ,ACT (4) ,AKP(4) ,ACP{4) 

COM MON/AEEBP/B J (4) 

COM MON/IFT/ITP 

COM MO N/I S M 3/ MS 

COM MO N/I R/I H 3 4 

CCM MO N/C YMl /CY5 6 , C Y5 7 

COMMON/NGAT/N WFC ,NP6 , NTRT 

COM MO N/LIN A/II NA , LAFI 

COM MON/CYM 2/C Y 58 

CCM MO N/NG AD/NHNI ,NCVF, NPPS 
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n n a 


C 

c 


COM MON/VMAKE/V MAXI 

DATA BLANK, MINUS/1 H ,1H-/ 

NC-TJT = 6 
NSIZRY = 225 

BEAD (5,902) NP3 ,NE4 , NP5,NP6,NTET 

RE AD ( 5 , 902) NB , NSP, MAXN , NES , NF, M ABR G, NAS , IGYRO, JY FC 
WRITD (6,910) 

IF ( JYEO, 3Q, 0) GO TO 9 
IF (MAXN. IF, 1 ) GO TO 9 
MAXN-1 
WRITE (6,9*17) 

9 CONTINUE 

IF (MP3. EQ« 0) GO TO 1 3 

WRITS (5, 9 1 2)b»1, NN2 , N N3 , NN4 , NN 5 , NN6 , N WF C , NH W I , N CV F , NPSS 
WRITS (6 , 91 3) I WAKE, NAERO, MXIT, KTKST , IH34 
WRITS (6,914) MP3 , NP 4 , NP 5 , N?6 , NTRT 

WRITS (6,916) NB , NS P, UAZN , NES , NF , MAERO, N A S , IG YRO , 0 YEO 

13 READ (5, 992) HNFAZ , NNBS, HFF ,IINA ,LAFI 
IF (NP3, E0, 0 ) GO TO 14 

WRIT'S (6,919) NNFAZ,NNB3,»FF,IINA,LAFI 

14 NFP1=NF+1 
MXSMI=2*NF+1 
MXKQ=12*MXSMI 

READ (5, 902) MFL AP,NFEA ,NCT ,NFLA P , NS , NBS, NPS , N3T , NIG 
IF (NP3. BO, 0) GO TO 15 

WRITE (6,917) MFLAP, N FE A , NCT , NFLAP, NS,N3S, NPS ,NST , NIG 

15 HCOLS = 6 
NCSB = 1 

READ (5, 902) MINPN , MA XPN 
NP = MINPN 

IF(IWAKB.EQ.C) GO TO 19 
IF (HAEiO.GT.C) GO TO 19 
WRITS (6, 9C 8) 

NP= MAXPN+1 
GO TO 210 
19 CONTINUE 

IF(NP3.F0.f ) GO TO 16 
WRITE (6,911 ) MINER’ , MAXPN 


COMPUTE ARRAY INDEXING VARIABLES 


16 NR 3D = 8+MFLAP 

MXT2P1 = (2*MAXN+1) *NSP 
MXCSB = 12*NCSB 
MXCPK = 12* NCOLS 
MXCPL = MXSMI*MXSMI 
NRIFC = MXT2P9 + N 3* NRBD 
MXCPK = MXC PL* MXCPM 


20 



c 

c 

c 


MXSHB = MXCPL* MXCS B 
UFA SB = ffXCPL* NCSE 
MXFA3 = MXCPL*NCGLS 
NB1FC = HRIFC* KXT2P1 
NEI SC = NEIFC*NBBD 
HEBC = MXCPM*MXSMI 


NE5BC=MXCS B*KXSMI 
MXQ = NRIFC*MXSMJ 
RXTK.N = NEIFC + NB*NBISC 

El AD (5, SOI ) CY4,C¥9,CYlI3,CYl4,CYl5,CY11,CY46,GY4 7,CY48,CY49,CY5Q, 
1CY54,CY55,CY56,CY58, VMAXI 
CY57=COS (CY56) 

CY56 = SIN{CY56) 

IF (8:P3« EO. C ) GO TO 17 
WRITE (6, 921) 

WRITE (6 ,918) CY4 ,CY9,CYl 3 , CY 1 4 , C Y15 , CY 1 1 , CY 4 6 , CY47 ,C Y48 , C Y49 , C Y50 
1CY54,CY55,CY56,CY58, VMAXI 


t 


READ IN SWASH-PLATE VARIABLES 


17 IF ( NSP. ECU 0) GO TO 30 
IF (JYEOvEQ. 1 ) GO TO 20 

READ (5, 901) S W G J , S WE I ,SWM,SWR 

IF (NP4. EQ.C ) GO TO 18 

WRITE ( 6 , 920) S WG J , S WEI , S WM , S W R 

18 CALL SUPPHA 

READ (5, 901) (A K ( J) ,J=1,NES) 

READ (5,901) (AC ( >7) ,J=1,NES) 

READ (5, 901) CAPK ,CAPC 

READ (5, 901) (AKT (J) , J=1 ,NES) 

READ (5,901) (ACT (J) , J=1,NES) 

READ (5, 901) (AKP(J), J = 1 ,N1S) 

READ (5,901) (ACP(J), J=1 ,NES) 

READ (5, 901 ) (BJ(J), J=1 ,NB3) 

GO TO 30 

20 READ (5, 901) GYM ,GYK, GYC , GIX ,GIY ,GIZ ,G KX , GKY ,GCX ,GCY, GRP 
IF (NP4, E0» 0 ) GO TC 22 

WRITE (6,922) GYM ,GIX ,G YK ,GIY ,GYC,GIZ,GKX,GKY, GCX , GCY , GRP 
22 G KB = , 5* (GKX+GKY) 

GKP = . 5*<GKX-GKY) 

GCB = . 5* (GCX+GCY) 

GCP = . 5* (GCX-GCY) 

30 IF (MAXN.GT.NF) GO TO 32 
ITP = 2*NF 
GO TO 34 
32 ITP=2*MAXN 
34 CALL PHASE 

READ (5, 901) (D MS (MS) , MS=1,NB) 

READ (5,901) (AKCI(MS), MS=1,NB) 

READ (5, 901) (TAU(MS), MS=1,NB) 
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non 


HEAD (5 , 9 )1 ) (S MLA (MS) , MS=1,NB) 

CY 2 = (G. , 1 . ) 

CY7 = CY4*CY4 
CY8 = 2.0*CY4*IGYBO 
CY12 = 3.1415927 
CYl 3 = CY) 2*. 25 
CYl 6 = SIN (CYl 5) 

CYl 7 = COS (CYl 5) 

CYl 8 = NAS 
CY40 = NS 
CY41 = NS 
CY42 = NFIA 
CY43 = NF 
CY44 = MXSKI 
CY45 = NP 
CY51 = COS (CY4 6) 

CY52 = SIN (CY46) 

READ IN AND COMPUTE BLADE STATION DATA 

DO 40 1 = 1 , NSIZSY 
40 Y (I) =0. 

IF (NPS.EO.fi) GO TO 11 
WRITE (6,955) (1, 1 = 1, 7) 

11 DC 2S 9 MS=1 , NB 
DO 120 J = 1 , N S 

READ (5,909) Ml , M2 , K3 , M4 , M5, MS , M7 
Y (201) =M 1 

Y (202) = M2 

Y (203) = M 3 

Y (204) =M4 

Y (205) =M 5 

Y (206) = M6 

Y (2C 7) =M 7 

IF(MI.EQ.O) GO TO 10 
CALL YA ( J) 

IF (NP5.EQ.C) GO TO 12 
1=1 
11 = 5 

WRITE (MOOT, 956) J, Ml , M2 , M3 , M4 , MS , M6 , M7, MINUS , I , Y ( I) , 

1 (BLANK , X, Y (2) ,1=2,4) , MINUS , 1 1 , Y (IT ) 

WRITE (HOUT, 957) (SINUS, I., Y (I) ,1 = 6,7) , (BLANK, I, Y (I) ,1=8, 10) 
WRITE (NOUT, 957) (BLANK ,1 , Y (I) ,1=75,79) 

WRITE (UOUT, 957) (BLANK, I, Y (I) ,1=18^,182,2) 

12 IF ( M2 . EQ . 0) GO TO 3 
CALL BYB(J) 

IF (NP5.EQ.fi 1 ) GO TO 3 

WRITE (NOUT, 9 57) (BLANK, I, Y { I) ,1 = 8,10) 

3 IF (M4.UO.O) GO TO 5 
CALL EYE ( J) 
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IF (NP5. EQ. 0) GO TO 5 
IF (Ml .NF.O. OR. M2.NE.0) GOTO 4010 
1=150 
11=181 

WHITE (MOOT, 9 56) J,Ml,M2, M3,M4 ,M5 , M6 , M7 , MI NUS , I ,Y (I) , 

1 (BLANK , I, Y (I) ,1=151,154) 

WHITE (NOUT , 957 ) BLAN K, II , Y (II) 

GO TO 5 
4010 1=150 
11=181 

WHITE (NOUT, 957) MINUS, I, Y (I) , (BLANK, I, Y (I) ,1 = 151,154) 

WRITS (NOUT, 957) BLAN K, II , Y ( II) 

5 IF (M5, EQ.O) GO TO 7 
CALL YEIGID 
WRITS (6,915) 

IF (NP5. EQ. 0 ) GO TO 7 

IF{Ml.NB.0.OR.U2.*ng.0.OR.M4.SE.0) GO TO 6020 

WHIT'S (NOUT, 9 56) J,Ml , M2 , M3 , M4 , M5 , M6 , M7, (MINUS, I, Y (I) ,1 = 184,186) , 

1 (BLANK, I, Y (I) ,1 = 187,188) 

11=189 

WRITE (NOUT, 957) BLANK, II , Y (II) 

GO TO 7 

6020 WHITE (NOUT, 957) (MINUS , I, Y (I) ,1 = 184,186) , (BLANK, I, Y (I) ,1=187,188) 
11=189 

WRITE (NOUT, 957) BLANK, II, Y (II) 

7 IF (M6. EQ.O) GO TO 90 
CALL YSK 

IF (NP5. EQ.O) GO TO 90 

IF(M1. NE,0,OE,M2.NE.O. OR.M4.NE. 0.OH.M5. NE.Q) GO TO 8000 
11=196 

WRITE (NOUT ,956 ) J,M1 , M2, M3 /M 4 , M5 , M6 , M7, (MINUS, I, Y (I) ,1 = 1 90 ,192) , 

1 (BLANK, I, Y (I) ,1 = 1 93,1 94) 

1 = 195 

WHITE (NOUT, 957) BLANK, I, Y (I) , MINUS, II, Y (II) 

GO TO 90 
8000 11=196 

WRITE (NOUT, 957) (MINUS, I, Y (I) ,1=190,192) , (BLANK, I, Y (I) ,1=193,194) 
1=1 95 

WRIT'S (NOUT, 9 57) BLANK, I, Y (I) , MI HUS ,11, Y (II) 

GO TO 90 

10 IF (M2. EQ.O) GO TO 3 
CALL BYB(J) 

IF (NP5 . EQ. 0) GO TO 3 
IF (MI.NE.P) GO TO 2005 

2000 WRITS (NOUT, 9 56) 0, Ml , M2 , 83 , M4 , H5 ,36 , M7 , (BLANK, 1, Y (I) ,1=8,1’)) 

GO TO 3 

2005 WRITE (NOUT, 957) (BLANK, I , Y (I) ,t = 8 , 10 ) 

GO TO 3 

90 DO 200 L=1 , NSI ZE Y 
M- (J-1) *NSIZEY + L 
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200 SD(H)=Y(L) 

DO 203 1=1 , N SI Z>i Y 
203 Y (I) =0. 

12ft CONTINUE 

IF (MAERO . FQ. 0) GO TO 206 
DO 119 K=1 , NAS 
DG 119 L=1 , NFF 

CPSY (K,I.) =COS (L* ( (K-1) * (2,*CY12) /NAS) ) 

119 5PSY (K,L) =SI N ( L* ( (K-1) * (2. *CY12) /NAS) ) 

CALL A j5 BOM 

208 CONTINUE 

STORE BLADE STATION EAT A ON DISK 

WRITS (2) SD 

209 CONTINUE 
REWIND 2 
REWIND 3 

210 CONTINUE 
RETURN 

901 FORMAT (8F10. f) 

902 FORMAT (1515) 

907 FORMAT (3X,*HAXN MUST BE .LI. 1 IF JYRO IS 1, M AX N IS RESET TO 1*') 

908 FORMAT (3X,*IF IWAKE = 1, THEN MAERO MUST 33 INPUTTED AS 1*) 

909 FORMAT (811) 

910 FORMAT (THl , 50. X , 32HNASA HELICOPTER DYNAMIC RESPONSE/) 

911 FORMAT (32X,8H MINPN =,13,119 MAXPN =,I3) 

912 FORMAT (32X, 811 NNl = ,I4,10H NN2 =,I4,1QH NN3 = ,I3, 

1 11 H NN4 =,I3,11H NN5 =,I3/32X,8H NN6 =,I3, 

2 11 H NWFC =, 13 , 1 1 H NHWI =,13,11 H NCVF =,I3, 

3 11 H NPBS =, 13) 

913 FORMAT (32X, 8H IWAKE =,I3,1'1H NAS BO =,I3,11H EXIT =,I3, 

1 1 1 H KTBST =,I3 , 1 1 H IH34 =,I3) 

914 FORMAT (32X, 8H NP 3 =,I3,11H NP4 = ,l3,11 T i N?5 = ,!?, 

1 11 H N P6 =, 13 , 1 1 H NTRT =,I3) 

915 FORMAT (3.X, *IF MACRO NOT = 0 RIGID SECTIONS MUST BE AVOIDED, MAERO 
1= 0 CENTRIFUGAL FORCE Y (187) MOST BE INPUTTED, */ 3 X,* SUGGEST USE OF 

2 STIFF ELASTIC SECTIONS INSTEAD*/) 

916 FORMAT (32X, 8H NB =,I3,11H NSP =,13,918 HAXN =,I3, 

1 11 H NES = , I 3 , 1 1 H NF = , 13 , /3 2X, 3TI MAERO =,I3, 

2 11 H NAS =, 13,118 I GYRO =,I3,11H JYRO =,I3) 

917 FORMAT (32X,8H MFLAP =,13,118 KFEA =,I3,^1H NOT =,I3, 

1 1 1 H NFL A P = , 1 3 , 1 1 H NS =,I3,/32X,8H N3S =,I3, 

2 11H NFS = , 1 3 , 1 1 H NET =,13,118 NIG =,I3) 

918 FORMAT (3 (/) , 29 X, 1 7HOPER ATI NG SPEED =,'F10.4 , 19X, 

1 1 6H SPEED OF SOUND = , Fl 3 . 4/3 3X, 1 6HDPSSIT Y OF AIR =, 

2 FlO» 4 ,1 1X,2 4 H AIRCRAFT FORWARD SPIED =, FlO. 4/28X, 

2 18HROTOR TILT ANGLE =, F10 . 4 ,17X,1 3HGBA VIT Y CONSTANT = , F1C. 4/43X, 

3 6H BETA =, Fl C»4, 21 X,14H ALPHA CYCLIC =,F10. 4/35X , 1 1HROLL RATE =, 
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4 F10.4,17X # 18HCOLLBCTIV3 ANGLE = , Fl 0 . 4/3 2X , 1 4HHQTQR RADIUS = , 

5 Fl 0 • 4 , 2 1 X , 1 4 H PULL UP RATE = , Fl C . 4/32X, 1 4HT UBNI NG RATS = ,F10. 4, 
617X ,184SINfi (BANK ANGLE) = , F 10 . 4/3 3X , 1 3HR00T LENGTH =,F10.4, 
719X/16HVHL (I, J) LIMIT = ,F10.4) 

919 FORMAT (32X,8H NNFAZ =,T3,11H NN3S =,I3,11H NFF =,13, 

1 11 H IINA =,I3,11ff LAFI = ,I3) 

929 FORMAT (2 (/) ,31 X, 1 5 HS WAS H PLATS GJ =, F10. 2 , 23 X , 1 5HS WftSHPL ATE SI =, 

1 F10, 2/29X, 17HSWASHPLAT3 MASS = ,F10 » 4 ,1 6 X , 1 9 HSN A3 HP LATE' RADIUS =, 

2 FIT. 4) 

921 FORMAT (50X, 27HHELICOPTEH ROTOR PARAMETERS//) 

922 FORMAT (2 (/) ,35X,11HGYS0 MASS =, Fl 0 . 4 , 26 X , 9HGY30 IX = , Fl 0. 4/1 9 X, 

1 27 HCOLLECTIVS GYRO STIFFNESS = , Fl 0 • 2 , 2 6 X , 9 H GYP 0 I Y =,F10. 4/21X, 

2 25HCOLLSCTIVE G YE 0 DAMPING = , F 10 , 4 , 2 6X , 9HG YRO .12 = , ?19 . 4 , /22 X , 

3 24HLATERAL GYRO STIFFNESS = , F10 . 2 , 1 4X, 2 1 HLONG GYRO STIFFNESS =, 

4 F10# 2/24X, 22HLATSRAL GYP 0 DAMPING =,F10.4,16X, 

5 19HL0NG GYRO DAMPING = , F1U 4/3 3X, 13HGYR0 RADIUS =,?10. 4) 

955 FORMAT (3 (/) ,8H SECTION, 7(2H M , 1 1 ) , 3 7 X , 2 1 H I N PUT Y VECTOR VALUES) 

956 FOR M AT (1 5 ,3X ,713 , 5 (1 X , Al , 2H Y ( , 13 , 3H) =,E10,3)) 

957 FORMA T (29X , 5 (1 X, Al , 2HY ( , 13 , 3H) =,E10.3)) 

END 
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SUBROUTINE SDPPHA 
RIAL CHI (4) 

CO MMO N/RNAMEl/CS 1 (4,6) ,SN1 (4,6) 

COMMON- /NO 3 /NS , NSIZSY , NFEA , NF.S, MAXN , NFPl 

READ (5, 700) (CHI(J), J=1 ,NES) 

IF(MAXN.EQ.O) GO TO 13 

MXNT2=2*MAXN 

DO 12 J=1 , NES 

DO 12 NL=1 , MXNT2 

AEG =NL*CHI (J) 

CSl (J,NL) =COS (AEG) 

12 SHI (J,KL) = SIN (AEG) 

13 RETURN 

700 FORMAT (8 FI 0,6) 

END 
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SUBROUTINE PHASE 
REAL PHIM (4) ,PHIMP (4) 

COMMON/RNAME/CS (4,20) ,SN (4,20) 

COM MON/I FT /I TP 

COM M0N/N04/N COLS , NB , N E, NEIFC,.NSISC, MXTKN , NI 

READ (5,700) (PHIM (MS) , MS = 1 , N3) 

READ (5,700) (PHIMP(MS) , MS = 1 ,NB) 

IF (ITP.EQ.O) GO TO 11 
DO 10 MS = 1 ,NB 
DO 10 L=1 ,20 

ARG = L* (PHIM ( MS) - PH IMP (MS) ) 

CS (MS ,L) = CCS (ARG) 

10 SN ( MS , L) = SIN (ARG) 

11 RETURN 

700 FORMAT (8F10.6) 
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SUBROUTINE Yh (IS) 


Y A- 

MA£S MATRIX FOR BLADE- ELEMENTS OF VECTOR -Y- TO EE CREATED. 

COM MO N/CYR/CY4 , CY7 ,073, CY9 , CY10 , CY 1 1 , CY 1 2 , C Y 1 3, CY 1 4 , C Yl 5,C Y1 6 , 

% CYl 7, CY19 , CY20 ,CY21 ,C Y22 ,C Y23 ,CY24 ,CY25 ,CY26 ,C Y27 , CY28 ,CY?9, 

2 CY30 ,CY31 ,CY32,CY33,CY34,CY35, CY36 , CY37 ,CY 38 , C Y3 9 ,C Y4 6 ,C Y47 , C Y4S , 

3 CY49 , CY50, C Y51 , CY52 ,€ Y53 

COMMON /NO 3 /N S , N SI ZE Y , N FEA , NS S , MA XN , NF Pi 
COMaSN/SUB/Y (225) 

READ (5 ,7C 1 ) (Y (I) ,1 = 1 ,10) , (Y (K) ,K=75,79) ,Y (180) ,Y (182) 

IF (IS-NFEA) 590,60 7,600 
590 Y(10) =Y{10) +CY49 
600 Y (1 ) = - Y (1 ) 

Y<5) = - Y (5) 

Y<6) = -Y(6) 

Y (7) = -Y (7) 

Y (1 1 ) = SIN {¥(3)) 

Y (12) = COS ( Y ( 8) ) 

Y (1 3) = Y (1 1 ) * Y <1 1 ) 

Y (1 4) = Y (1 2) * Y (12) 

Y (15) = SI N (Y (9) ) 

Y (16) = COS (Y (9) ) 

Y (1 7) = Yfi5) * Y (1 5) 

Y (1 8) = Y(16) * Y { 1 6 ) 

Y (1 9) = SIN (Y (10 ) ) 

Y (20) = COS (Y (10) ) 

Y (2 1 ) = Y (1 9 ) * Y (19) 

Y (22) = Y (20) * Y (20) 

Y (23) = Y (1 1 ) * Y (1 2) 

Y (24) = Y (11) * Y (1 5) 

Y (25) = Y (1 1 ) * Y (16) 

Y (26) = Y (1 1) * Y (1 9) 

Y (27) = Y (1 1 ) * Y ( 20 ) 

Y (2 8) = Y ( 1 2) * Y <1 5) 

Y (29) = Y (12) * Y (1 6 ) 

Y (30) = Y (1 2) * Y ( 1 9) 

Y (3 1 ) = Y (1 2) * Y (20) 

Y (32) = Y ( 1 5) * Y (1 6) 

Y (3 3) = Y (1 5 ) * Y (19) 

Y (3 4) = Y (1 5) * Y (20) 

Y (3 5) = Y (1 6 ) * Y (1 9) 

Y (3 6) = Y (16) * Y (20) 

Y (37) = Y (1 9 ) * Y (21) 

Y ( 3 8) = Y (32) * Y (20) 

Y (3 9) = Y (32) * Y (1 9) 

Y (40) = Y(4) * Y (1 1 ) + Y (3) * Y (12) 

Y (4 1 ) = Y (3) * Y (1 1 ) - Y (4) * Y (12) 
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Y (42) = Y (18) * (Y (22) -Y (21) ) 

Y (43) = Y (18) * Y (37 ) 

Y (44) = Y (18) * Y (21 ) + Y (1 7) 

Y (45) = Y (17) - Y (18) * Y (21) 

Y (46) = Y ( 1 8) * Y (22) +Y(17) 

Y (4 7) = Y (1 8) * Y (22) - Y (17) 

Y (48) = Y (1) * Y (2) 

Y (4 9) = Y(1) * Y (2) * Y (2) 

Y (50 ) = Y (5 ) + Y (4 9) 

Y (51) = Y (5) - Y (6) 

Y (52) = Y (5 ) - Y (7) 

Y (53) = Y (6) - Y (7) 

Y (54) = Y (51 ) + Y (4 9) 

Y (55) = -Y (53) + Y (49) 

Y (56) = Y(5) - Y (53) + (2. * Y(49)) 

Y (57) = - Y (52) - Y (6) 

Y (58) = Y ( 7) - Y (51) 

Y (80) = 0. 

5KTUPN 

701 FOFM&T (8F10.6) 



SUBROUTINE BYB(IS) 

C BEND MATEIX- 

C 

COMMON/CYR/CY4,CY7,CY8,CY9,CY10,CY11,CY12,CY13,CY14,CY15,CY16, 

1 CY17,CY19»CY20,CY21 # CY22,CY23, CY24 , CY25 , CY26 , CY27 ,€ Y28 , C Y29 , 

2 CY30,CY31,CY32,CY33,CY34,CY35,CY36,CY37 ,C Y38 , C Y39 , CY46 , CY47 , CY48 , 

3 CY49 ,CY5C # CY51 , CY52 / CY5 3 

COfl HOD /NO 3/NS , NSI ZEY ,NFKA , NBS , MAXN , NFPl 
COMMGN/SAIN/SD (742 5) 

COMMON/SUB/Y (225) 
r 

R1AD(5,912) Y ( 8) , Y (9) , Y (10) 

IF (Y (?S-2) .SQ.2. ) GO TO 600 
IF ( IS— NFE5A) 5 9 0 , 6 C: € , 6 0 0 
590 Y (10) = Y (10) +CY49 
6( 0 IF (Y (2.01) -1.) 70 0, 735,700 
700 Y (11) = SIN (Y (8) ) 

Y (12) = COS ( Y ( 8) ) 

Y (15) = SIN (Y (9) ) 

Y (1 6) = COS (Y (9) ) 

Y (1 8) = Y (16) * Y (16) 

Y (1 9) = SIN (Y (10) ) 

Y (20) = COS (Y (10) ) 

CM.MM1 1=ABS (Y (11) ) 

CJJH Ml 2= ABB (Y (12) ) 

CMF Ml 5=ABS ( Y (1 5) ) 

CBN Ml 6 = A BS (Y (16) ) 

CMMMl 9=ABS ( Y (19) ) 

CMMM20-ABS (Y (2v) ) 

IF(CMMM11 .LT. .IB-04) Y(11)=D.O 
IF (CMMMl 2 „ LT. . IB-34) Y (1 2) -0, <> 

I F ( CP MM 1 5 .LT. .11-04) Y (15) =0.0 

I F (CF MMl 6 . LT. . IE-04) Y (16) =0.0 

IF ( CM MM 1 9 .LT. . IT-9 4) Y (19) =0.0 

IF (CMMM2S . LT. . IB-04) Y (2 ) =0, 3 


Y (21) 

-= 

Y ( 1 9 ) 


Y ( 1 9 ) 

Y (22) 


Y (20) 

4; 

Y ( 2C; ) 

Y (2 3) 

= 

Y (1 1 ) 

# 

Y (12) 

Y (24) 

- 

Y (1 1) 


Y (15) 

Y (25) 

■= 

y r i ) 


Y (16) 

Y (2 6) 

- 

Y( * 1) 

¥ 

Y (19) 

Y (27) 

■cr 

y r i ) 

¥ 

Y (20) 

Y (2 8) 

= 

Y (1 2) 


Y (15) 

Y (2 9) 


Y (12) 

4s 

Y(16) 

Y (3 9) 

— 

Y { i 2) 

•?< 

Y (1 9) 

Y (31) 

- 

Y CL) 

if- 

Y (26) 

Y (33) 


Y (i 5) 

m 

Y(19) 

Y (34) 

=r 

Y (1 5) 


Y ( 30) 

Y (3 5) 

= 

Y (16) 

'jft 

Y ( ' 9 ) 

Y (3 6) 

= 

Y (3 6) 

4F 

Y(2v) 
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735 1= (IS-2) *NSIZEY+11 
NSP 1=N5+T 

IF (IS.EQ.1) GO TO 741 
IF(IS.EQ.HSPI) GO TO 741 

739 Y (60) = SD(I) 

I = 1+1 

Y (61) = SD (I) 

I = 1 + 3 

Y (62) = SD (I) 

I = 1+4 

Y (63) = SD (I) 

I = 1+6 

Y (64) = SD (I) 

K = I + 6 
1=1+1 

M = 6 5 

DO 740 L = I,K 
Y(M) = SD (I) 

740 M = M+1 
1=1 + 7 

Y (71 ) = SD (I ) 

1 = 1+1 

Y (72) = SD (I) 

1 = 1 + 1 

Y (73) = SD (I ) 

1 = 1+1 

Y (74) = SD (I) 

GO TO 742 

741 Y (61 ) =1. 

Y (68) =1. 

Y (70) = 1. 

Y (74) = 1. 

742 BET DEM 

912 FORMAT (3 FI 0.6) 

END 
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SDBBOIJTIWE BYE (IS) 

COM MO H/C YR /C Y 4 , C Y7 ,C Y 8 ,C Y9 , C Yl 3 ,C Yl 1 , C Yl 2 , C Yl 3 , CY 1 ft , C Y1 5 , CYl 6 , 

1 C Yl7,CYl9 # CY?i,CY2 1 , CY2 2, CY23 , CY24 ,CY2 5, CY26, C Y27 ,C Y28 ,C Y29 , 

2 CY39,CY31,CY32.,CY?3,CY34 ,CY35,CY36 ,CY37 ,CY38 , CY3 9 , CY46 ,0147 , CY48 , 

3 CY49,CY5'‘?,CY51,CY52,CY53 

COMMON /N03/NS,NSIZEY ,NFEA , HNS, MAXN,SFPl 
CdMHOS/S AIN/S L (7425) 

COMMUVSUR/Y (225) 

HEAD (5,701) < Y (I) ,1 = 1 53,154) ,Y (181) 

IP (Y (2^4) -1 . ) 650,648,6 51 
648 1 (150) =■ - Y (150) 

650 IF ( Y (204 ) - 1 . ) 9 ,9 , 1 1 
9 Y (1 54) = C. 

IF (IS , BO. 1) GO TO 11 
ISM 1 =IS“ 1 
no 13 J2=1 , ISMl 
J 2 M 1 = (J2-1) ♦NSIZKY 

10 Y (154) =Y (154)-SD(J2M1+1) *SD (J2M1+3) *CY7 

11 Y (1 55) =-Y{150) * S 0 S T (Y (154) / Y (157) ) 

Y (1 56) = - Y (1 50) * SQRT (Y (154) / Y (i 53 ) ) 

Y (1 57) = SINH (Y (155) ) 

Y (153) = COSH (Y (155) ) 

Y (159) = SINH (Y (156)) 

Y (1 6 3) = COSH (Y (156)) 

Y (161) = Y (155) ^ Y (155) 

Y (162) = Y (161) * Y (161) 

Y (163) = Y (156) * Y (156) 

Y (164) = Y (163) * Y (163) 

IF (Y (15 5) -.Cl) 1570, 3 570 ,15 51 
1550 Y (16 5) = Y (157) / Y(155) 

Y (166) = (Y (157) -Y (155) ) /{Y (161) *Y (155) ) 

Y (167) = Y (158) 

Y (168) = (Y (158) - 1. ) / Y (161) 

IF (Y (156) -.01) 1 580,1533,1 560 

VY.GT..1 VZ.GT.1 

1560 Y (169) = Y (159) / Y(156) 

Y (1 70) = ( Y (159) — Y (1 56) ) / ( Y (1 56) * Y(163)) 

Y (171) = Y (160) 

Y (172) = (Y (160) - 1.) / Y (163) 

GO TO 652 

GAMMA Y )= .1 

1570 Y (1 65) = 1. + Y (161) / 6, + Y (1 62) / 120. 

Y (166) = .166666 + Y(161) / 12I.+Y(162) / 5043. 

Y (1 67) =1. + Y (1 61) / 2. + Y (1 62) / 24. 

Y (163) = .5 + Y (161) / 24. + Y(162) / 720. 
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IF (Y (155) —.01) 1580,1585,156-) 

1580 Y (1 69) - 1. + Y (1 63) / 6. + Y <1 64) / 12 0. 

Y <170) = . 166666 + Y (163) / 120. + Y (164) / 5540, 

Y (171) = 1. + Y (163) / 2. + Y (164) / 24. 

Y (1 72) = .5 + Y (163) / 24. + Y(164) / 720. 

GO TO 652 

651 Y (158) =1, 

Y (1 60 ) =1 , 

Y (165) =1. 

Y (166) =.16 666666 

Y (167) =1 . 

Y (168) =.5 

Y (169) =1. 

Y (170) =. 1 6 66 666 

Y (171) =1. 

Y (172) = .5 

652 RSTUSN 

701 FORMAT (8F10.6) 

END 
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SUBROUTINE YRIGID 

THE Y VALUES - -18 4-1 8 6— ARE DELTA X 

COHHOH/SUB/Y {225) 

C 

READ (5/701) (Y (I) ,1=1 84,189) 

RETURN 

701 FORMAT (8F10.6) 

END 


Y, ANDZ 
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SUBROUTINE YSK 


-Y- VECTOR FOE STIFFNESS (SK) MATRIX 


COMMON/SUB/Y (225) 

C 

READ (5,701) (Y (L) , L=190 ,1 96) 
C TO INVERT 

IF (Y (206) -1 • ) 79,78,79 

78 Y (190) = -Y (190) 

Y (191) = -Y (191) 

Y (192) = -Y (192) 

Y (196) = -Y (196) 

79 RETURN 

701 FORMAT (8F10.6) 

END 



CVFRLAY (NS6,1 ,2) 

PROGRAM LIFT 

REAL LIF (10, 24) 

REAL LTH(lfs) 

REAL LTKMP 

COMMON /C YE/ CY4,CY7 ,C Y8 ,C Y9 ,CY1G ,C Yl 1 , C Yl 2 , CYl 3 , CYl 4 , CYl 5 , CY 1 6 , 

1 CY17,CY19,CY20,CY21,CY22,CY23,CY24,CY25,CY26,CY27,CY28 ,C-Y29, 

2 CY30,CY3?,CY32,CY33,CY34,CY35,CY36,CY37 ,C Y38 ,CY39,CY46 ,CY47 , CY48 , 

3 CY49,CY50 ,CY51 ,CY52,CY53 
COM MO N/IAR 1 /1AR (10) 

COMMON/I WA1 /IWAKE, NAEHO, MA350 

COM MO N /NO 5/ N P, N SP , MODI , MFP , NOUT , M XQ , N AS , NBS , NET , 1* PS 
COMMON / W Al / GTBST - 

COKMON/WA2/GAM (240) 

HEAL GAMMA (10,24) 

COMMON/WA3/VZA1 (13,33) ,VYAl (13,33) ,VZPl (13 ,33) , VYPl (13,33) 
CCMMON/WA4/ALPH1 (10,24) ,VELF (13,24) , LTH 
COM MO N /I W A 4 /N N 4 , NN 5 , N N 6 
COM MON /NG AD /NH WI ,NCVF , NPR S 

IF (IAR (1) ,BQ. 1) GO TO 5 
DO 6 Kl =1 , NAS 
ANGA = 2,0*<Kl-1)*CYl2/NAS 
C Y2 0 = SI N ( A N G A) 

CY21 =CQS (ANGA) 

CY53=-CY47* (CY21*C Y51 +C Y2G*C Y52) 

CCY53 =C0S (CY53) 

SCY53=3IN ( C Y 5 3 ) 

DO 6 11 =1 , NAEFO 
12 = IAR (II) 

122 = 12-1 

CALP = COS(ALPHl (I1,Kl) ) 

SALP=SI N ( AL PHI (1 1 , Kl) ) 

CALPA=CCY53*CALP+SCY53*3ALP 
S AL PA=CCY53*SALP-SCY53*CALP 
LIE (II, Kl) = (VZAl (1 ,12) -VZAl (1,122) ) *CALPA + 

1 (V YAl (1 ,12) - V Y A 1 (1,122) ) *SALPA 
IF(NH'WI.LE.O) GO TO 6 
DO 6 Jl =1 , NHWI 
<32 = Jl +1 

LIF (II , Kl) =LIF (II, Kl) +2.0* ((VZAl (J2,I2) *COS ( ANG A* Jl +VZ Pi ( 32 , 12 ) ) 

1 -VZAl ( J2,I22) *CGS (ANGA*Jl +VZP1 (J2, 122) ) )*CALPA+ 

2 (VY Al (J2,I2) *COS (ANGA* Jl +VYP1 (J2,I2) ) 

3 -V YAl ( 32 , 122 ) *COS (ANGA*Jl+ V YPl (J2,I22) ) )*5ALPA) 

6 CONTINUE 

IF (NN4 , EQ. 0) GO TO 30 

WRITS (6,103) ( (LIF (I , J) , 1=1, NAESO) ,J = 1,NAS) 

WRITE (6, 104) ( (VELF (I , J) ,I=1,NABRO) ,J=1,NAS) 
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30 GDI =0 * 

G-d. 

DO 21 15=1 , NAEEO 
DO 21 J5=1 , NAS 
K5= (J5-1 )*NAEB0+I5 

21 GAMMA (15, J5) =GAM (K5) 

DO 3 13 = 1 , NAEP.O 

DO 3 K3=1,NAS 

X5 = CY10* VSLF (13, K3) *LTH (13) *CY4*CY5O*CY50 
IF (X5» NS. 0.0) GO TO U 
LTERP=0» 0 
GO TO 9 

4 LTEMP=LIF(I3,K3) /X5 

9 GDIFF=LTEHP-GAHHA(I3,K3) 

GDI=GBI+ABS ( GDIFF) 

GAMMA (13 , K3) =GABMA (13, K3) +. 01*NCVF*GDIFF 
G=G+A8S (LTEMP) 

3 CONTINUE 

I F ( G . GT . • f ( C Of 1 ) GO TO 10 
« RITE (6,905) 

GTE ST =0.0 
GO TO 12 
10 GTEST=GDJ/G 
12 IF (MN5. FQ.'D ) GO TO 31 

WE IT 33 (6 , 102) ( (GAMMA (I, J) ,1=1 ,NAERO) , J=1 ,NAS) 

102 FOF MAT (* POUND CIRCULATIONS */(8Gl4.5)) 

103 FORMAT (* LIFT ON BLADE */(8Gl4. 5)) 

104 FOE MAT (* FREE STREAM VELOCITIES */(8Gl4.5)) 

31 DO 22 16=1 , NAEEO 
DO 2.2 J 6 = 1, NAS 
K6= (,76-1 ) *NAERO+I6 

22 GAM (K6) =GAMMA (16, J6) 

RETURN 

5 WRITE (6,101) 

101 FORMAT (104HC AERO MATRIX NOT PERMITTED BY CONVENTION AT BLADE TIP S 
IT AT ION, REDEFINE INPUT AND REPLACE BY MASS MATRIX) 

905 FORMAT (3 X, ^CIRCULATIONS ARE LESS THAN . IE-05, GTE ST SET = S.C*/) 
STOP 
END 
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OVERLAY (NS6,1,3) 

PROGRAM WAKE 
C 

INTEGER H 
C 

DIMENSION SIG (240) 

C 

COMHON/VEL1/VEL (10 ,24) 

COMMON/WA2/G AM (240) 

COB BO H/IS A 1 /I W A KB , N ABRO f BABB 0 
COMMGN/T WA2/ITI, NANR, IMAX, KTEST 
CO M. MO N/I W A 4/NN 4 , NN 5 , S N 6 
COB MON/KG AT/NWFC ,NP6 , NTRT 
COMMON /IWA3/ MXIT 
COM MO N/VMAKE/V MAXI 
C 

C IF FIRST TIME, READ GAMMAS OR SET TO CONSTANT 

C 

MSB 1=0 

NONE=0 

N2=2* NAERO 

NPl =NAERO+1 

N AS=NANR/NAERO 

IF (ITI.GT.0) GO TO 3 

IF (KTEST. EQ. 0 ) GO TO 1 

IF (ITI.BO.O) READ (8) (GA B (K) , K=1 , NANR) 

GO TO 3 

1 IF(NWFC.EQ.O) GO TO 18 

READ (5,901) (GAK (K) , K=1 , NANR) 

GO TO 3 

18 DO 2 K=1 , NA NS 
2 GAM (K) =.02 
C 

C LOOP OVER TOTAL NUMBER OF AERODYNAMIC POINTS 

C 

3 DO 8 K = 1 , N.AN8 

READ (9) MS2T, (SIG (M) ,M=1, NANR) 

GO TO 6 
C 

C COMPUTE INDEX VARIABLES FOB MATRIX MULT 

C 

6 J= ( MSBT-1) /NALRO+1 
I=J*NABRO-*-1--MSET 
C 

C DO MATRIX MULT, WHERE H ORDERS THE GAMMAS 

c 

SUM=0. 

DO 7 M = 1 , N A S E 

H=N2* ( (M-D/MBRO) +NP1-M 

IF (NWFC.BQ. 1) GO TO 7 


38 



non 


IF (ITI.EQ.0) H=M 

7 S0M = SIG (H) *GAM(H) +SUH 
N0NB=1 

SET I,J INDUCED VELOCITY FOB NEXT ITERATION 

VEL (I ,J) *-SDH 

IF ( VM AXI . EQ . 0 , 0 ) GO TO- 8 

V ABS= ABS (VEL (I, J) ) 

IF ( VABS . GT. VMAXI) VEL (I , J) = VBL(I,J) *VMAXI/VABS 

8 CONTINUE 

IF (NN4 .20. 0 ) GO TO 21 
WRITE (6,108) (GA'H(N) , N=1 ,NANR) 

10 8 FORMAT (* GAMMA FROM WAKE #/{8Gl4.5)) 

21 IF (ITI. NE, MXIT) GO TO 19 
IF {NP6.E0.0) GO TO 19 
WRITE (4,901) (GAM (K) ,K = 1 , NANS) 

19 CONTINUE 

WRITE (6, 101). ( (VEL (I , J) ,1=1 , NABEO) , J=1 , NAS) 

101 FORMAT (* INDUCED VELOCITIES FROM WAKE */(8Gl«.5)) 
901 FORMAT (5G16. 7) 

IF (KTBST.NE.C?) REWIND 8 

REWIND 9 

RETURN 

END 
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OVERLAY (WS6,1 , 4) 

PROGRAM F A h RAY 

INTEGER CY3 # CY18,CY«C, CY41, CY 42 , CY4 3, CY 4 4 , C Y4 5, CY1 8A 
I NT EGER FAN (4) ,JSM(6) 

HEAL ERE AL (144) 

REAL VZ A 1 (13,33) ,VYA1 (13,33) ,VZP1 (13,33) ,VY?1 (13,33) 

COMPLEX AM A (1700) 

COMPLEX IMF (1230) 

COMPLEX TH (7 2) ,FP (12) ,R (144) ,T (72) ,D (144) ,C (144) 

COMPLEX B (1 800) ,H(60) ,SML3{300) ,SMLC(3C0) ,SHLD.(3C0) 

COMPLEX FAB (15 3 ) , FLB (150) ,CTB (1 50) 

COMPLEX FAI! ( 5) ,FLH(5) ,CTH (5) 

COMPLEX FAS B (25) ,CTSB(25) ,FLSB{25) , FLSC (25) ,FASD(25) ,CTSD(25) 
COMPLEX BSAVB (180'3) ,HSAVB(63) , SMLBS V (300 ) ,SMLCS V ( 3C0) 

COMPLEX SMLDSV (30C) ,SHAPE(12) 

COMPLEX EPS (123) ,DETSV 
COMPLEX CY 1 , C Y 2 , C Y 5 , C Y6 
COMPLEX FASC (25) 

COMPLEX EXPOS (5) 

COMMON/AMFT/AMF 
COM MO N/ A MAT/ A MA 
COM MON/S OB/ Y (2 25) 

COM MCN/J3 PS A/EPS , LETS V 

COM MO N/BTSl /B , SMLB , S MIC , SMLD 

COMMON/BTS2/FAE, FLB, CTB 

CO MMON/BT S 3/FASB ,C1 SB , PL SB , F LSC , FA SD , CT SD , F AS C 
COMMON/RET S/H , FAH, FLH , CTH 
COM MON/IP Hi /NN FAZ , NNB S 
COM MON/I SI /NFF 

COM MON/C YI/CY3,CY1 8 , C Y41 ,C Y4 1 , C Y42 , C Y43 , C Y4 4 , CY45 ,CY18A 
CCMMOW/C YC/CYl ,CY2,CY5,CY6 

COM MON/CYR/CY4 ,C Y7,C Y8,CY9 ,CY10 ,CYl 1,CY12,CY13,CY14,CY1 5-,CYl6, 

1 CY17,CY19,CY20,CY21,CY22,CY23,CY24,CY25,CY26,CY27,CY28,CY29, 

2 CY30,CY31,CY32,CY33,CY-34,CY35,CY36,CY37 ,CY38 , CY39 , CY46 ,CY47 , CY48 , 

3 CY49 , CY 50 , C Y 5 1 , C Y 5 2 , C Y 5 3 
COM MON/I SMB /MS 

COM MON/NOl/MXCPK , KXKQ , MXSMB, MFASB, MXFAB 

CCMMCN/N02/MXCSB ,NLSBC , MXCPM,NEBC 

COM MO N /N03/NS , NSI ZEY , NFL A, N1S , MAXN , NFPl 

C0MM0N/N04/NC0LS.,NB ,NF, NBIFC , NEISC , MXTKN ,NIG 

CCM M0N/N05/N P, NSP, MODE, HFP, NOUT , MXQ , NAS , N BS , NET , NPS 

C0BMCN/NQ6/NRBD, NR IFC ,M XSMI , NFLAP, MFLAP , NCT , MXT2P1 , NCS B 

CCMM0N/N08/J YEO 

COMMON/WA1/CONVG 

COMMON/WA5/ALLI 

C0MM0N/INA2/ITI,NANF.,IMAX,KTSST 
COM MQN/IRA3/MXIT 
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COM MON/SSI N/SD (7425) 

C 

COMMON/I WA1 /I WIRE, NAEB.O, MAEfiO 
COMMON/WA3/VZAl,VYAl ,VZPl, VYPl 
COM MON/ENAME/CS (4 ,20) ,SM (4,20) 
COMMON/NGAD/EHVI ,KCVF , MPES 
C 

c 

DATA EH/1, 15, 29, 42, 57, 70/ 

DATA FAM/4 , 8 , 1 2/ 

C 

M X N F F=N F F* 2 + 1 
NSD=NS*N5IZEY 
BEAD (?) SD 
NC€T= 1 
NCFBA=1 
NCFLP=t 
11=0 

5 DO 20 L = 1 , M XKQ 
20 H (I) =CMPLX <G. ,0 , ) 

DO 15 L= 1 , MXSMB 
SHLB (L) =CMPLX (0. ,0. ) 

SMLC(L) =CMPLX (0. ,0.) 

15 SMLD (L) =€MPLX (0. ,0 . ) 

DO 22 K=1 , MXCPK 
22 B(K) = (0.,C.) 

DO 30 1=1 , MXSHI 
L= (I“1 ) * MXS MI + 1 
K*L*MXCSB-9 
SMLC (K) = (1. ,0. ) 

K=K +1 

SHLB (K) = (1. ,0. ) 

K=K +6 

SMLD (K) = (1 , ,0. ) *MFLAP ' 
l Ml a (L-1) *MXCPM 
DO 30 M = 1 , 6 
K=LMl +EM (H) 

30 B (K) = (1 * ,0 • ) 

DO 817 1=1 , MFASB 
K= (1-1) *12+3 
817 FASC (I) =S MLC (K) 

IF (MFP-.EQ.O) GO TO 14 
DO 13 JX=1, MXSHI 
LX=JX-NP-NFPl 
LXA=I ABS (LX) 

IF (LX) 11,10,12 

10 EXPOS (JX) = (1. ,0. ) 

1 ' GO TO 13 

11 U1=CS (MS,LXA) 

U2=SN (MS , LXA) 
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EXPOS (JX) =CMPLX (01 ,-U2) 

GO TO 13 

12 U 1 = C S ( M S , L X ) 

U2=SN (MS, LX) 

EXPOS (JX) =CMPL X (U 1 ,U2) 

13 CONTINUE 

14 CONTINUE 

DO 501 IS=1, NS 

ISM 1 = (IS-1) *NSIZEY 

NCPLX=1 

DO 50 L =1 , NSIZIiY 
M=I SMI + L 

50 Y (L) =SD (M) 

IF (Y (205) ) 80,80,55 
55 CALL RIGID (RBBAL) 

LGO=1 

60 DO 79 1=1 , MXSMI 
IMl = (1-1) *12 

51 DO 61 M = 1 , 1 2 
K=I Ml + M 

61 TH { M) =H ( K) 

GO TO (62,63) , NCPLX 

62 CALL MLRC1 (RKEAL,TH) 

GO TO 64 

63 CALL MLCCl (R , TH) 

64 DO 65 M=1 , 1 2 
K = I M 1 +M 

6 5 H (K) =TH (M) 

52 DO 79 J=1, MXSMI 
L= (J-1 ) * MXSMI + 1 
LMl = (L- 1 ) *MXCPM 
DO 66 M=1 , MX CP M 
K=LMl +M 

66 T ( M) = B ( K) 

GO TO (67,6 8) , NCPLX 

67 CALL MLBC2 {RREAL,T, NCCLS) 
GO TO 69 

68 CALL MLCC2 (R,?,NCOLS) 

69 DO 70 M=1 , MXCPM 
K=LM1 + M 

B (K) -T(M) 

70 CONTINUE 

GO TO (77,71) , NCFX A 

71 LBMl = (L-1 ) * MXCSS 
DO 72 M= 1 , M XCSB 
K=LBMl + M 

72 TH ( M) =SMLC ( K ) 

GO TO (73,74) , NCPLX 

73 CALL HLEC2 (RRKAL ,TH,NCSB) 
GO TO 75 
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74 CALL MLCC2 (R,TH,NCSB) 

75 DO 76 M=1,MXCSB 
K=LBM1+M 

76 SMLC(K)=TH{M) 

77 GO TO (78,711) ,NCCT 

711 LBMl = (L-1) * ; MXCSB 
DO 721 K=1,MXCSB 
K=LBM1+M 

721 TH (M) -SMLB (K) 

GO TO (731,741) ,HCPLX 

731 CALL MLRC2 (RREAL ,TH, NCS B) 
GO TO 751 

741 CALL NLCC2 (R,7H,NCSB) 

751 DO 761 M~1 , MXCSB 
K=LBMl + M 

761 SH.LB (K) =TH (M) 

78 GO TO (79,712) ,NCFLP 

712 LBW1= (L-1) *MXCSB 
DO 722 M=1 , MXCSB 
K=LBMl+M 

722 TH (M) =SHLD (K) 

GO TO (732, 742) ,NCPLX 

732 CALL HLRC2 (RREAL, TH, HCSB) 
GO TO 752 

742 CALL MLCC2 (B,TH,NCS3) 

752 DO 762 MM, MXCSB 
K=L B Hi + M 

762 SMLD(K) =TfI (M) 

79 COST! NOB 

GO TO (80, IOC, 120) , LGO 
.80 IP (Y (202) ) ICC ,100,85 
85 CALL BIND (RREAL) 

LGO = 2 
GO TO 60 

100 IF (Y (204 ) ) 1 20 , 120 ,102 
102 CALL BLAST (H) 

LGO =3 
NCPLX = 2 
GO TO 60 

120 DO 200 1 = 1 , MX Si'll 
I Ml* (1—1 ) *1 2 
KSML=I-NFPl 
CY3=KSML 

CY5*CY1-CY.2»CY3*CY.U 
CY6=CY5*CY5 
KM.NP=KSML-NP 
IF (Y (206) ) 145,145,125 
125 CALL STIFF (R,CY5) 

LG 0=1 

130 DO 132 M =1 ,12 



K=in +h 
132 TH (M) =H (K) 

CALL MLCC^ (K ,TH) 

DO 135 M=1,12 
K = I M 1 + M 

135 H(K)=TH{M) 

127 DO 4 04 

L= (J-1) *3XSMI + I 
LMl = {L-1) *MXCPM 
DO 136 M=1,MXCFK 
K=LMl +M 

136 T(H)=B(K) 

CALL MLCC2 (B.,T,NCOLS) 
DO 137 M = 1 , MXC PM 
K=LMl +M 
B (K) =T (M) 

137 CONTINUF 

GO TO (142 ,1 38) , NCFEA 

138 LBMT= (1-1) *MXCS8 
DO 139 M=1 , MXCSB 
K=LBMl +M 

139 TH (H) = SMLC (K) 

CALL KLCC2 (B ,TH,NCSB) 
DO 140 M=1 , MXCSB 
K=LBMl +M 

140 SMLC (K) =TH (M) 

142, CONTINUE 

GO TO {400 , 4(1) , NCFLP 

401 LBMl= (L-1) * .MXCSB 
DO 402 M=1, MXCSB 
K =L B M, 1 + M 

402 TH (fl) =SflLD (K) 

CALL MLCC2 (R , TH , NCSB) 
DO 403 M=1 , MXCSB 
K=LBMl+M 

403 SMI.D (K) =TH (M) 

400 CONTINUE 

GO TO (404,40 5) ,NCCT 
415 LBH1 = (L-1)*MXCSB 
DO 4{ 6 M=1, MXCSB 
K=L BMl +M 

406 TH (M) =SMLB <K) 

CALL MLCC2 . (B.,TH,NCSE) 
DO 407 M=1 , MXCSB 
K=LBMl +M 

407 SMI, 3 ( K) =TH (M) 

404 CONTINUE 

144 GO TO {145,171), LGO 

145 I F { Y (201) ) 200,200,150 
150 CYl 8 A=- KMNP 


44 



165 CALL MASSB(L,FP) 

L GO =2 
GO TO 130 
171 DO 180 K -i ,12 
K-IMl +H 

180 H (K) =H (K) +FP (M) 

2C0 CONTINUE 

IF {Y (203)) 204 ,305 , 204 
204 IF(IL.GT.O) GO TO 2X9 
IF (I NIKE . EQ. f?) GO TO 
HEAD (3) AM F 
GO TO 209 

2,01 PE AD ( 3) AMA , AMP 

209 II=IL+1 

GO TO (2141,21 0 ) , N Z FEA 

210 DO 214 K=1 , MX SMB 
SMLCSV(K) -S MLC (K) 

214 SMLC (K) = (0. ,<5* ) 

2141 GO TO (2142,211) ,NCFLP 

211 DO 212 K=1 , MXS MB 
SMLD5V (K) =SMLB (K) 

212 SKID (K) = (0. ,0. ) 

2142 GO TO (2143, 281) ,NCCT 

281 DO 282 K=1 , MX.S MB 
SMLBSV(K) =S MIF (K) 

282 SML3 (K) = (0 » ,fl . ) 

2143 DO 215 K=1 , MXCPK 
BSAVE (K) =B (K) 

215 B (K) = (C.,C.) 

DO 216 K=1 , MXKQ 
HSAV® (K) =H (K) 

216 H (K) = <0, ,0, ) 

NT'IMS=M.X.SM.I + 1 

DC 3( 0 N Q = 1 , N F P 1 
CYl 8A=NQ-1 
NTIM3 = NTIMS-1 

nshft=nq-nfpi-1 

NBACK=1 

217 CALL A.EROB (C , D , FP , IL) 

DO 250 J 0*1, NT I MS 
NSHFT=NSHFT+1 

C Y3=NSHFT 

219 IQ*JQ+NQ-NBACK 
I=JQ+ NBACK-1 
IQMl = (XQ.-1) *12 
I Ml = (1-1 ) *12 . 
CY5=CYl-CY2*CY4*NSHFT 
DO 220 M=1 , 1 44 

220 H (M) =CY5*C (M) +D (M) 

DO 22 5 M=1 , 12 



K=IQM1+M 

2,25 TH (M) =HSAVK (K) 

CALL MLCC1 (R ,TH) 

DO 228 H=1 ,12 
K=IM1+M 

228 H (K) =TH (M) +H (K) 

IF ( NS HFT . N L . 0) GO TO 235 
KK= (IL-1)*MXNFF*3 + (NP + NFPl - X+NFF) *3 
DO 230 M=1 , 3 
K= (1-1) *12+FAM (M) 

230 H(K) =H(K) +AMF(KK+H) 

235 DO 250 J=1 , MXSMI 

L= (J-1) *MXSM I+I 
LQ= (J-1) * MXSMI + IQ 
LH1 = (L-1) *BXCFM 
LQWl = (LQ-1 ) +MXCPM 
DO 240 H=1,HXCPM 
K=LQM1+M 

240 T(M) =BSAVE(K) 

C AIL MLCC2 f B , T , SCOLS ) 

CO 245 M=1,MXCPM 
K = L Ml + M 

245 B (K) = T(M) + B (K) 

GO TO (249,246), NCFEA 

246 L BQ = (LQ-1) *HXCSB 
1B = (L-1) * flXCS 13 
DO 247 M=1,MXCSB 
K=LBQ+M 

247 TH (M) =SMLCSV <K) 

CALL MLCC7 (E ,TH,UCSB) 

DO 24 8 M=1 , MXCS B 
K=LB+M 

248 SSLC (K) =TH (N) +5HLC (K) 

249 GO TO (231,261) ,SCFLP 

261 LBQ=(LQ-1)*MXCSB 
L 8= (L-1) *MXCS8 
DO 262 f*=1,MXCS3 
K=LBQ+M 

2 62 TH (M) =Si1LDSV (K) 

CALL M LC C2 ( R , T H , ?F CS B ) 

DO 26 3 M =1 , M XC SB 
K=LB+H 

26 3 SMLD(K) =T?T ( M) + SMLD (K) 

231 GO TO { 2 5t ,232), SC CT 

232 LBQ = (LQ-1) XCSB 
LB= (L-1) *MXC5B 
DO 233 H =1 , MXC3 B 
K=LB0+M - 

233 TH (H) =SBL3SV (K) 

CALL MLCC2 (R ,TH, SCSB) 
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DO 234 M=1 , MXC SB 
K=LB+M 

234 SML9 ( K) =TH ( M) + SMLB (K) 

250 CONTI NOE 

251 CONTINUE 

IF(NQ-I) 300,300,255 

255 GO TO (256,300) ,N3ACK 

256 CY18A=-CY18A 
NBACK=2 
NSHFT=-NFPi 
GO TO 217 

300 CONTINUE 
305 CONTINUE 

IF (MFE . 30. 1) GO TO 350 
IF (IS-NFEA) 326,310,326 
310 NCFEA-2 

DO 320 1=1 , MXFAB 
K= (1-1) *12 + 4 

320 FAB (I) =B (K) 

DO 321 1=1 , MXS MI 
K = (1-1) *12 + 4 

321 FAH (I ) =H (K) 

312 IF (IS-NCT) 323,222,322 

322 DO 324 I=1,HFASB 
K= (1-1) *12 + 4 

324 FASB (I) =SMIB (K) 

323 IF (IS-NFLAP) 326,325,325 

325 DO 327 1=1, KFASB 
K= (1-1) *12 + 4 

327 FAS D (I) =SP!I.D (K) 

326 CONTINUE 

IF { I S- N CT ) 34 (.:■ ,326,340 

328 HCCT=2 

DO 329 1=1, MXFAB 
K = (1-1) *12 + 3 

329 CTB (I ) =B (K) 

DO 330 1=1 , NXSPil 
K= (1-1) *12 + 3 

330 CTH (I) =H (K) 

DO 331 I=1,KFASB 
K = (1-1) *12 + 4 

331 CIS B (I) = SMI.P (K) 

• IF (IS-NFLAP) 323,332,332 

332 DO 334 I=1,MFASB 
K = (1-1) *12+3 

334 CTSD (I) = SMI,D (K) 

333 CONTINUE 

340 IF (IS-NFLAP) 500,335,506 

335 NCFLP=2 

DO 336 1=1, MXFAB 



K= (1-1) *12+11 

336 FIB (I) = B (K) 

DC 337 1=1 , MX SMI 
K = (1-1) *12 + 1 1 

337 FLH (I ) =H (K) 

316 IF(IS-HCT) 339,338,338 

338 DO 342 1=1 , MFASB 
K= ( I- 1 ) *12+1 ^ 

342 FI, SB (I) -SFLB (K) 

339 IF (XS-NFKA) 506,343,343 

343 DO 344 1=1 ,MFASB 
K= (1-1) *12 + 11 

34 4 FLSC (I) =SMLC (K) 

GO TO 5M 

350 IF (IS . I!0* NFEA) NCFEA-2 
IF (IS • IiQ. NFLAP) NCFLP=2 
IF (IS.SQ. NCI) HCCT = 2 

IB = 1 

IF (NCT ,GT. IS) 1 0=0 
IC=1 

IF (NFS A . GT, IS) IC=0 
ID=1 

IF (NFLAP ,GT. IS) ID=0 
DO 351 IR=1 ,12 

351 SHAPE (IS) =C MPLX (3. ,0 . ) 

J3-NF*72 

J2=RF*12 

<J 1 = M X T 2 P 1 + (MS-1) * HERD 

357 DO 366 J=1,MXSMI 
JMl =J-1 

I Jl =JM1*TJ5?BC + J3 
I J2=JMl*NBSBC +J2 
I Kl = J Ml * S 31 FC+ Jl 
DO 366 IB =1,1 2 
K=I Jl+IR 
L=I J2+IR 

SHAPE (IR)=SHAPK(IR)-IC*SHLC(L) *EPS (IKl+7) *BXP08 (J) -IB*SMLB (L) * 
1BXPOB (J) *EFS (I.Kl+8) 

IF (ID .FQ. 0) GO TO 367 

SHAPE (IE) =SHAPB (IB) -ID*SHLD (I) *£PS (IK1+9) *EXPOB (J) 

367 DO 366 IP=1,6 
I1=IK1+IP 
K1=K+(IP-1) *12 

366 SHAPE (IR) =SHAP3(IR) +B (Kl) *EPS (Ll)*BXPOB (J) 

DO 36 8 15=1 ,12 
L3=J2+IR 

358 SHAPE (IR) = S H APE (IR) +H (L3) 

IF ( IWAKE . BQ . 0) GO TO 4 
IF { MS . GT. 1 ) GO TO 4 
IF(NP.GT.O) GO TO 2 
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VYA1 (1,IS) =-b.ML (SHAPt (8) ) 

V ZA 1 (1,TS) =RLAL {SHAPE (1 2) ) 

VYP1 (1,1$) =0.0 
VZPl (1,IS) =0.0 
GO TO 4 
2 MNPl=NP+i 

V YAl (BHP1,IS)=-CABS (SHAP3C8) ) 

VZAl (MNPl , IS) =CABS (SHAPE (12) ) 

Xl = SB AL (SHAPE (8) ) 

X2=AIMAG (SHAPE (8) ) 

X3=REAL (SHAPE (12) ) 

X4=AIMAG (SHAPE (12) ) 

VYP1 (MNPl ,IS) =0.0 
VZPl (MNP 1 , IS) =3.0 

IF(X1.EQ.0, .AND.X2.EQ.0.) GO TO 3 

V YPl (MNPl, IS) = ATAN 2 { X2 , Xl ) 

.3 IF (X3 . EQ. 0 . . AND. X4 . IQ. 0 . ) GO T 0 4 
VZPl (MNP 1 , IS) =ATAH2 (X4,X3) 

4 CONTI N HE 

499 CONTINUE 

IF (I WAKE . FQ. O) GO TO 603 
IF(NPHS. N.S, 6) GO TO 603 
IF (ITI-MXIT) 500 ,603,60 3 
603 IF (MS, GT. T) GO TO 808 
IF (IS. GT. 1 ) GO TO 80 9 
KMl=NF*NBJFC+1 
KM2 = MK1 + MXT2 Pi -1 
IF ( NS P . EQ. 0) GO TO 80 3 
WRITS (6,950) 

IF ( JYEO .EQ. 1) GO TO 804 
IF ( MAXN-1 ) ■ 8 C 1 ,802,303 

801 WRITE (6, 951 ) (EPS (IBP) , IBP-KMl , MM2) 

GO TO 808 

802 WRITE (6, 952) (BPS (IBP) ,IBP=fSMl , MM?) 

GO TO 808 

803 WRITS (6, 95 3) (EPS (IBP) , IEP = K Ml , MM2) 

GO TO 3 Cl 8 

804 IF ( MAXN . EQ* 1 ) GO TO 806 

805 WRITE (6,954) (BPS (ISP) , I2P=MMl , MM2) 

GO TO 80 8 

806 WRITS (6,955) (BPS (IBP) , ISP=MMl ,F M2) 

898 IF(IS.GT.I) GO TO 809 

WRITS (6,961) MS 
809 WRITE (6,991) IS 

WRITE (6,919) (SHAPE (IR) ,IR=1,12) 

500 CONTINUE 

501 CONTINUE 
RETURN 

901 FORMAT (//59X, 11 HSECT ION I = ,I3) 

950 FORMAT (/55X , 22HSWASHPLAT2 DEFLECTIONS/) 



951 F0RMAT(/6X,*W (0) =*, 2 (1PE12 . 4) /) 

952 FORMAT (/5X,*W (-1) =*, 2 { 1 PEl 2, 4) , 6X, *« (0) =* , 2 (1 PE12. 4) , 5X, 

=*,2 (1PB12.4)/) 

953 FORMAT (/5X , *MAXN .GT. 1, W (I) WHERE I GOES FROM -MAXN TO +MAXN*/ 
2 (2X , 10PE1 3.4) ) 

954 FORMAT {/6X , *Z (0) =*,2(1 PEI 2.4)/) 

955 FORMAT(/5X,*PHIXR = *, 2 (1 PEl 2. 4) ,6X,*Z (0) =*, 2 <1 PE12. 4) , . 
15X,*PHIYR =*,2 (1PE12. 4) /) 

961 FORMAT (//51 X , 29HSTA7E VECTORS ON BLADE NUMBER, 13) 

910 FORMAT (/1 6 X,*UX =* , 2 ( 1 PEI 2 . 4) , 6X , *N =*, 2 (1PE12. 4) , 3X, 

1 *PHI X =*,2<1PE12.4)/17X,*T =*, 2 (1 PEI 2. 4) ,5X,*U Y =*, 

22<1PE12.4) ,3X,*PHI Z =* , 2 (1 PE12. 4) /I 6X,*MZ =* , 2 (1 PEl 2. 4) , 

34 X VY =* f 2 (1 PE 12. 4) , 5X,*--JZ =*,2 (1PF12.4) /13X, 

4*PHI Y =*, 2 (1PE12. 4) ,5X,*MY =* , 2 (1 PEl 2 . 4 ) , 6 X , 

5*VZ =*,?. (1PS12-.4) /) 

END 
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non 


SUBROUTINE BEND (B) 

Y (60) THRU Y (74) = PREVIOUS STATION DATA ARRAY VALUES 
REAL B (144) 

DIMENSION SA (3 , 3) , SAT (3,3) , SB (3,3) 

COMMON /SUB/ Y (225) 

SA (1,1) = Y (29 ) 

SA ( 1 , 2) = Y (25) 

SA (1 ,3) = -Y (15) 

S A (2, 1) = Y (1 9) * Y(28) — Y ( 2 7 ) 

SA (2 ,2) = Y (33) * Y ( 1 1 ) + Y(31) 

SA (2, 3) = Y ( 35) 

S A (3 , 1 ) = Y (34 ) * Y (1 2) + Y(26) 

SA ( 3, 2) = Y (34) * Y (1 1 ) -Y(3l>) 

SA { 3 , 3 ) = Y (36) 

SAT (1,1) = Y ( 6 8) 

SAT (1,2) = Y (63) * Y (67) - Y (66) 

SAT (1,3) = Y ( 7 2) * Y (61) + Y(65) 

SAT (2,1) = Y (6 4) 

SAT (2,2) = Y (71) * Y (60) + Y(70) 

SAT (2,3) = Y (72) * Y (60) - Y{69) 

SAT (3,1) = - Y ( 52) 

SAT (3,2) = Y (7 3) 

SAT (3,3) = Y (7 4) 

DO 20 1=1,3 
DO 20 J = 1 , 3 
SB (I, J) =0.0 
DO 20 K = 1 , 3 

20 SB (I, J) = SB (I, J) +SA (I , K) *SAT (K,,J) 

DO 5 1=1,144 
5 B (I ) =0. 

11 3 (1 ) = SB (1,1) 

B (5) = SB (1,2) 

B (S) =-SB ( 1 , 3) 

B (1 4 ) = SB (1,1) 

B (20) =- SB (1,2) 

B (24) = SB (1 ,3) 

8(27) = SB (1,1) 

B(3C) = SB (1,3) 

8(34) = SB (1 , 2) 

B (40) = SB (1,1) 

B (43) = SB (1,3) 

B (47) = SB (1 ,2) 

B (49) = SB (2,1) 

B (53) = SB (2 ,2) 

B (57) =— S3 (2,3) 

3(63) = SB (3,1) 

E (66) = SB (3,3) 
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B (70) = SB ( 3 , 2 ) 
B{76) = SB (3,1) 
B (79) = SB (3,3) 
B (83) = SB ( 3 , 2 ) 
B (86) =- SB ( 2 , 1 ) 
B (92) = SB (2,2) 
B (96) =— SB ( 2,3) 

B (97) =-SB (3,1) 
B (101) =-SB { 3 , 2) 
B (105) = SB (3,3) 
B (111) = SB (2,1) 
B (1 1 4 ) = SB (2,3) 
B (1 1 8) = SB (2,2) 
B (124)= SB (2,1) 
B (127) = SB (2,3) 
B (131)= SB (2,2) 
B (1 44) = SB (3,3) 
B (134)= SB (3,1 ) 
B (1 40) =-SB (3,2) 
SET URN 
END 
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SUBROUTINE BLAST (E) 

ELASTIC MATH IX FRO M -Y- VICTOR 

COMPLEX E f 1 44) 

COMMON /SOB/ Y (225) 

DO 5 1=1,144 
5 1(1) =(0.,0«) 

1590 DO 1.6.00 K=1 ,144,13 

1600 B (K) =1. 

E (28) = Y (150) / Y (151) 

. E (54) = 7(15^) * Y (165) 

B (55) = Y (150) * Y (150) / Y (152) * Y(168) 

E (56) = Y (150) * Y (150) * Y (151) / Y(152) * Y ( 1 5 6 ) 

E (66) = Y (1 67) 

E (67) = Y(15C) / Y (152) * Y (16 5) 

1(68) = Y < 1 5 il ) * Y (1 50) / Y (152) * Y (1 68) 

E (7 8) = Y (154) * Y (150) * Y (165) 

E (7 9) = Y (16 7) 

E (80) = Y (150) * Y (165) 

H (1 06) = Y (1 5C ) * Y (169) 

E (1C 7) = Y (150) * Y (150) / Y (153) * Y (1 72) 

E (108) = Y (150) * Y (150) * Y (1 50) / Y(153) * Y{170) 
E (118) = Y (171 ) 

E (119) = Y (1 50) / Y (153) * Y <1 69) 

E (1 20 ) = Y(15C) * Y (150) / Y(153) * Y(172) 

R (1 30) = Y (1 54) * Y (1 50) * Y(16 9) 

S (1 31 ) = Y (171 ) 

E (132) = Y.(15C) * Y (169) 

RETURN 

END 
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SUBROUTINE RIGID (R) 

REAL R (144) 

COMMON/SUB/Y (225) 

DO 5 L=1 , 1 44 

5 R (L) =0. 

IF (Y (20 5) - 1 ) 80,79,80 

79 Y (1 84) = -Y (184) 

Y (1 85) = -Y (185) 

Y (1 86) = -Y (186) 

80 DO 10 L= 1, 144,13 
10 B(I) = 1 . 

R (6) =-Y (185) 

R (TO) = Y (186) 

R (44) =-Y(186> 

R (4 8) =-Y (185) 

R (51.) =-Y (186) 

R (54 ) = Y (184) 

P (7 4) = Y ( 18 5) 

R (80) = Y (184) 

R (99) =-Y (1 85) 

R (106) = Y (1 8 4 ) 

R (1 22) =-Y (1 86) 

R (1 32 ) =Y (184) 

R (3 9) = + (Y(186)*Y(189)+Y(185)* Y(188)) 

R (4 2) =- Y (184) *Y (189) 

R (46) =- Y (184) * Y (188) 

R (7 5) =~Y (186) *Y (1 87) 

R (7 8) = + ( Y (1 8 5) * Y (18 8) + Y(184) * Y (187) ) 

R (82) =- Y (186) * Y (1 88) 

R (123) = -Y (185) * Y (187) 

R (126) =-Y (185) * y (189) 

R (1 30) = + ( Y (184) * Y (1 87) + Y (186) * Y (189) ) 

RETURN 

END 
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SUBROUTINE STI FF (S K , CY5) 

COMPLEX CY5 
COMPLEX SK (1 44) 

COMMON /SUB/ Y (225) 

DO 5 1=1,144 
5 SK (I) = (G. ,0.) 

IF (Y (196) ) 198,199,198 

198 SK (3 9) = Y (196) 

GO TO 200 

199 SK (28) = Y ( 1 90 ) / (1.+ Y (193) * CY5) 

SK (67 ) = Y (1 92) / (1. + Y (19 5) * CY5) 
SK (1 1 9) = Y (191) / (1. + Y (194) * CY5) 

200 DO 201 I = 1,144,13 

201 SK (I) = (1. ,0. ) 

RETURN 

END 



SUBROUTINE ABHOB (AC,AD,VD,IL) 

C CREATE AERO DYNAMICS MATRIX FOR BLADE 

COMPLEX AC, AD, YD 
COMPLEX ASF (12 30) 

COMPLEX AHA (1700) 

INTEGER CY4Q , CY4 1 , CY42 
COMPLEX CY1 ,CY2,CY5,CY6 
INTEGER CYl 8 A 

INTEGER CY43,CY45,CY3,CY18,CY44 
DIMENSION AC {144) , AD{144) , VD{12) 

COMMON /AM FT/ AMF 

COMMON /AMAT/AMA 

COMMON/ISl/NFF 

COM MO N/C YC/C Y1 , CY2 , C Y 5 , C Y 6 

COMMON/CYI/CY3 ,C Yl 8 , C Y4 0 ,C Y41 , C Y42 , CY43 , C Y4 4 ,CY45 ,CY1 8 A 

COM MON/C YR/CY4,CY7,CY8,CY9,CY 10 f CYl 1 , CY 1 2 , C Yl 3, CY 1 4 ,C Yl 5 ,C Y1 6, 

1 CYl 7, CY19, CY20,CY21 ,C Y2 2 ,C Y23 , CY24 ,CY2 5,CY26,CY27,CY28 ,CY29 , 

2 CY39 ,CY31 , CY32, CY33 , CY34 , CY35 , CY36 , CY37 ,CY 38 , CY3 9, CY46 ,CY47,CY48, 

3 CY49,CY50,CY51,CY52,CY53 
MXNFF=2*NFF+1 

DO 10 J = 1 , 1 4 4 
AC (J) = (0.0, 0.0) 

10 AD (0) = (0.0, 0.0) 

DO 12 J=1 ,12 

12 VD { J) = (0.0, 0.0) 

IPQ = 34 * (IL-1 ) *CY44 + (C Y43 + CY1 8 A) * 34 
JPQ = 3^ (IL- 1 ) *M XNFF + (NFF+CYl 8A) *3 
19 FORM AT ( 1H ,* AC,AD,VD FROM AHA*/) 

IIK=IPQ+1 
I MK -IPQ +34 

21 FORMAT (1 HO , 10 El 2. 6) 

AC (39) = AMA (IPQ+1) 

AC (41 ) = AMA (I PQ + 2 ) 

AC (45) = AMA (IPQ + 3) 

AC (87) = A MA (I PQ + 4 ) 

AC (89) = AMA (IPQ + 5) 

AC (93) = AMA (I PQ+6 ) 

AC (13 5) = A MA ( IPQ + 7) 

AC (137) = AMA (IPQ •*■8) 

AC (141) = A M A { IPQ + 9) 

IF (CY18A) 9,8,9 

C SET DIAGONALS =1., CREATE MATRIX D. 

8 DO 101 I = . 1,144,13 
101 AD (I) =(1.,0.) 

9 AD (18) = AMA (IPQ+10) 

AD { 22) = AMA (IPQ+11) 

AD (37) = AM A (I PQ+1 2) 

AC (39) = AMA (.IPQ + 1 3) 

AD (41 ) = AMA (I PQ + 1 4) 

AD (42) = AMA (IPQ + 1 5) 
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AD (45) = AMA (IPQ + 1 6) 
AD (46) - AMA (I PQ+1 7 ) 
AD (82) = AM A (I PQ+1 8) 
AD (85) = A MA (I PQ + 1 9) 
AD (87) = AM A (IPQ + 20) 
AD (89) = AMA (1 PQ+2 1 ) 
AD ( 90) = ASA (IPQ+22) 
AD (93 ) = AMA (IPQ + 23) 
AD (94) = AMA (IPQ + 24) 
AD (126 ) = AMA (IPQ+25) 
AD (133) * A MA (IPQ+26) 
AD (135) = AMA (IPQ+27) 
AD (137) = AMA (IPQ + 28) 
AD (138) = AMA (IPQ+29) 
AD (1 41) * AMA (IPQ+30) 
AD (142) = A MA (IPQ+31 ) 
CDS ATE VECTOR D. 

VD (4) = AMF (JPQ+1) 

VD (8) =AMF(JPQ + 2) 

V D ( 1 2 ) = A MF (JPQ+3) 

BSTURN 

END 
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SUBROUTINE MAS5B (A, F) 

C 

C MASS MATRIX *A* FOR BLADE 

C 

INTEGER CY18A 
INTEGER €Y4G,CY4i,CY42 
INTEGER CY4 3 , CY45 , CY3 , CYl 8 , CY44 
COMPLEX CYl,CY2,CY5 f CY6 
COMPLEX A (144) 

COMPLEX F (1 2) 

COMPLEX CN9 

C CREATE MASS MATRIX FROM -Y- VECTOR FOR BLADE 

COM MON/SUB/ Y (225) 

COMMON/CYI/CY3 ,CY1 8,CY40»CY41 ,CY42 » CY43 , CY44 ,CY45 ,CYl8A 
CCMMON/C YC/CYl -,CY2 ,CY5, CY6 

COMMON/C YR/CY4 ,€Y7,C Y8 ,CY9 ,C Y10 ,CYl 1 ,CYl 2 ,CYl 3, CY14 , CYl 5, CYl 6, 

1 CYl 7 ,CYl9,CY2tJ,CY21,CY2 2,CY23,CY24,CY25,CY26 / CY27,CY28,CY29, 

2 CY30,CY31,CY32,CY33,CY34,CY35,CY36,CY37 ,CY38 ,CY39, CY46.,CY47 , CY48 

3 CY 49 , CY5C> C Y5 1 , CY 52 , C Y 5 3 
■C0KMON/CYM1 /CY56,CY57 
COM MON/C YK/C Y 5 4 , CY 55 

COM MO N/CYM2/CY58 
CSAVE-CY4 

TUF.N=CY48*CYi7 + CY55*CY57*CYlS+CY4 

RCLL=CY48*C Yl 6-C Y55*CY57*C Yl 7 

PtfLL=CY'54+CY55*CY56 

CY4=TURN 

C Y7 = CY4*'CY4 

CY8-2.0*C¥4 

CNl = Y (1)*CYl1i*CYl6 

HME=Y (3) -CY58 

EPHY=Y (2) + Y (4) 

r* 

v 

CY35 = Y (28) *Y (19) -Y (27) 

CY36=Y (31) + Y (24}*Y (19) 

CY37 = Y (26) + Y (28) *Y (20) 

CY3 8 = Y (24) *Y(20)-¥ (39) 

DO 5 1=1,144 
5 A (I ) = (D . ,0. ) 

DO 11 1=1,144,13 
11 A (I ) = (1 . ,0 . ) 

A (13) = Y { 1 ) * (CY6 - CY7 * Y ( I S ) ) 

A (15) = Y (48) * (CY8 * Y(35) * CY5 - CY7 * Y (39) ) 

A (17) = -Y(1) * (CY8 * Y { 36) * CY 5 + CY 7 * Y(39)) 

A (18) = -Y (48) * (CY6-CY7 * Y(18)) 

A (21) = — Y (1 ) * (CY 0 * Y (35) * CY 5 - CY 7 * Y(38)) 

A (37) = -Y { 48 ) * (CY8 * Y(35) * CY5 + CY7 * Y(38)) 

A (39) = Y ( 50) *CY6 + Y{55) *CY7*Y (42) +Y (48) *CY7* (HUE* Y (39) 

1 + Y (4) *Y (46) +CY58*CY35) -CN1*Y (2) *Y (35) 

A (4 1 ) = ~Y (48) * (CY8 * Y(15) * CY5 - CY7 * Y(43)) 
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A (42) =Y (56) *CY8*Y (35) *CY5*. 5+Y (54) *CY7* Y (38) 

A (45) =-Y (48) * (CYb-CY7*Y (44) ) 

A (46) =Y(57)*CY8*Y{36) *CY5*. 5+Y ( 52) *CY7*Y(39) 

A (73) = -Y ( 48) * (CY6 - CY7 ★ Y(18)) 

A (75) =-Y (56) #C Y8*Y (35) *CY5*. 5+Y (55) *CY7*Y (38) 

A (77) = Y (48) * ( CY8 * Y(36) * CY5 + CY7 * Y (39) ) 

A (7 8) = ( Y (49) +Y (7) ) *CY6 + Y(54) *CY7* Y (45) + Y (49) *CY7* (HMT5* Y (39) 

1 + Y (4) *Y (46) +CYS8*CY35) -CN1*Y (2) *Y (35) 

A (81) = Y (48) * ( CY8 * Y(35) * CY5 - CY7 * Y(38)) 

A (8 2) =- Y (58) *C Y8*Y (1 5) *CY5*. 5+Y (52) *CY7*Y (43) 

A (8 5) = -Y (1 ) * (CY8 * Y (36) * CY5 - CY7 * Y { 39 ) ) 

A (87) = - Y (48) * ( CY8 * Y (1 5) * CY5 + CY7 * Y{43)) 

A (89) = -Y (1 ) * (CY6 - CY7 * Y (46) ) 

A (90) = Y (48) * (CY8 * Y (36) * CY5 - CY 7 * Y(39)) 

A (93) = Y (1 ) * ( CY8 * Y (1 5) * CY5 + CY7 * Y(43)) 

A (1 23) =-Y (57) *CY8*Y (36) *CY5*. 5-Y (55) *CY7*Y ( 39) - Y (48) *CY7* (HME*Y (1 8 
1) +Y (4 ) *Y (39) +CY58*Y (29) ) -CN1*Y (2) *Y (15) 

A (1 26) = Y (58) *CY8*Y (1 5) *CY5*. 5 + Y (54) *C Y7*Y (43) -Y (4 8) *CY7* (HMB* Y (38) 
1-Y (4) *Y (43) +CY58*CY37) +CNl*Y (2) *Y (36) 

A (139)= Y (6) * CY6 - Y(52) * CY7 * Y(47) 

A (133) = -Y (1 ) * (CY8 * Y (35) * CY 5 + CY 7 * Y(38)) 

A (135)= Y (43 ) * ( CY6 - CY7 * Y (44) ) 

A (137) =-Y (1) * (CY8*Y(15) *CY5-CY7*Y (4 3) ) 

A (138)= Y (48) * (CY8 * Y (35) * CY 5 + CY 7 * Y (38) ) 

A (141)= — Y ( 1 ) * (CY6 - ZY7 * Y(44)) 

DO 6 1=1,12 
6 F (I) = (0 « ,0 « ) 

CYl8A=-CYl8A 
IF (CY18A) 81,3,81 

81 IF (CY18A-1) 82,15,82 

82 IF (CY18A+1) 190,23,109 

8 F (2) =-Y(1) *CY7* (SPH Y* Y (39) +HMB* Y (1 8) +CY58*Y (29) ) 

1 -CN7*Y (15) 

F (4 ) = Y (55) *CY7*Y (43) +Y (48) *CY7* (Y (4) *Y (43) ~HMB*Y ( 38) -CY58*CY37) 

2 +Y (2) *CNl*Y (3 6) 

F (7) = Y (54) *CY'7*Y (39) +Y (48) *CY7* (Y (4 ) ^ Y (39) + HH3*Y (18) +CY58*Y (/ 9) ) 

7 +Y (2)*CN1*Y (15) 

F (8) = Y ( 1 ) *CY7* (EPHY* Y (46) + HHE*Y (39) +CY58*CY35) 

2 -CN1*Y (35) 

F <1 1 ) = -Y{52) * CY7 * Y (38) 

F (7 2) =Y (1) *CY7* (EPHY*Y (4 3) ~HMB*Y (33) - CY 58* I Y 3 7) 

2 +CN7*Y(36) 

GO TO 700 
20 CONTI NfJB 
S=7 . 

25 CN2 = Y (1) *CY11*CY17*-. 5 

CN3 = ( CY 3 6*BOLL+C Y35* PULL) *10 SH 
CK4=(CY35*BOLL-CY36*PULL) *TUSN 
CN5= (Y (25) *EOLL + Y (29) *PULL) *TURN 
CN6 = ( Y ( 29) *EOLL-Y (25) *P(JLL) *TURN 



CN7= ('CY38*BGLL+C Y37*PULL) *Y?JRN 
CN8= (CY37*BOLL-CY38*PDLL) *TUBN 

CH9=Y (1)* (fiPHY»CS3+HHB*CN5+Cl58*PDLL*T0BH-S : *CY2*( 

1EPH Y*CN4 +H ^E*C'i36 + CY 5 8*HGLL*TUBN)-) 

F (2) =-CN9* Y (15)-CH2* (Y (25) -S*CY2*Y (29) ) 

F{8) =-CN9*Y (35) +CN2* (CY36-3*CY2*CY35) 

F (12) =+CN9*Y (36) -CN2* (CY38-S+CY2*CY37); 

F (4 ) = +Y (2) *F (12) -Y (53) *Y (36) * (CH3-S*CY2*CH4) +Y ( 57) *. 5* (-CH6 
1-S*CY2*CH5)r 

F (7) = -Y (2) *F (2) +Y (51) *Y (15) * (CM 3-S*CY2*CN 4) -Y (5 8) *.5* (-CN8 
1-S*CY2*CM7) 

F (1 1 ) =~Y (52) *Y (1 5) *(CN7-S*CY2*CH8) +Y (58)*, 5* (CN4+S*CY2* CM3) 

. GG TO ICO 
15 CONTINUE 

S=-1. 

SO TO 25 
ICO CY4=CSAVE 
CY7=CY4*CY4 
GY9=2.‘)*CY4 
RE TUB N 
END 
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SUBROUTINE MLR Cl (R,T) 
REAL E (144) 

COMPLEX T (72) ,THLD (12) 
DO 50 1=1,12 
50 THLD (I) = (0. ,0. ) 

DO 100 M =1 , 1 2 
( 8 - 1 ) *12 
DO ICO 1=1,12 
LM= MMl +1 

100 THLD (M) =R (LM) *T (I) +THLD 
DO 200 1=1,12 
200 T (I) =THLD (I) 

RETURN 

END 



SUBROUTINE HLEC2 (E ,T ,NCOLS) 
BEIL R (144) 

COMPLEX T (72) , THLD (72) 
NCR=1 2*NCOLS 
DO 50 1=1, NCR 
50 THLD (I) = (0. ,0. ) 

DO IPO N=1 , NCOLS 
NNl = ( N- 1 ) * 1 2 
DO ICO M=1 ,12 
MB1= <M-1)*12 

k=eni+h 

DO 100 1=1,12 
LM=MM 1 +1 
LN= NNl + 1 

100 THLD (K) =R(LM)'*T(LN) +THLD (K) 
DO 200 1=1, NCR 
200 T (I) =THI.D (I) 

SET URN 
END 
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SUBROUTINE MLCC1 (R,T) 

COMPLEX E (1 44) , T {72) ,THLD <1 2) 
DO 50 1=1,12 
50 THLD (I) = (0. ,0. ) 

DO 100 M =1 ,12 
M Ml = (H-1) *12 
DO 100 1=1,12 
LM= MMl +1 

100 THLD (M) =E(LM) *T(I) +THLD (M) 

DO 200 1=1,12 
200 T (I) =THLD (I) 

RETURN 

END 



SlJBFOUfl Hi’ MLCC2 (R,T,HCOLS) 
COMPLEX T (7 2) , THLD(72) # R{144) 
NCR =1 2 * N CO LS 
DO 50 1=1, NCR 
50 THLD (I) = (0. ,0, ) 

DO 10 fi N=1 , NCOLS 
HN1 = CM— 1 ) *12 
DO ICO Pi =1,12 
PJ Ml = (B-1>*12 
K =L N 1 +M 
DO 100 1=1,12 
LM=MM1+I 
LN=NNl+I 

100 THLD (K) =R <LM) *T (LM) + THLD (K) 

DO 200 1=1, NCK 

200 T (I) =THLD(I) 

RETURN 

END 
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OVERLAY (NS6 ,1,5) 

PROGRAM F K N S 
INTEGER P 
INTEGER ALP {6) 

INTEGER CY40 , C Y4 1 , C Y42 ,0 Y4 3 , CY4 4 ,C Y45 ,C Y3 , C Yl 8 , C Yl 8 A 

COMPLEX CYl ,CY2 ,CY5,CY6 

COMPLEX H ( 6 3 ) ,FAH{5) ,FLH(5) ,CTU (5) 

COMPLEX PKN (123) 

COMPLEX BXPOF 
COMPLEX ZKM 

COM MOM/C YI/CY3,CY1 8 , C Y4 0 , C Y4 1 , C Y42 , CY43 , C Y4 4 , CY45 , CYl 8 A 
CO M MO N/C YC/C Yl ,CY2,CY5, CY 6 

COM MO N/C YR /C Y4 , C Y7 , C Y 8 ,C Y9 , C Yl 4 , CYl 1 , C Yl 2 , C Yl 3 , CY 1 4 , CYl 5 , CYl 6 , 

1 CY17,CY19,CY20,CY21 , CY2 2, CY23, CY.24 , CY2 5 ,CY26, C Y27, C Y28 ,C Y29 , 

2 CY 30 , CY 31, CY 32, CY3 3, CY34,CY35,CY36,CY37 , C Y 3 8 , C Y3 9 , C Y4 6 , € Y 47 , CY 4 8 , 

3 C Y4 9 , C Y50 , C Y5 1 , CY 5 2 , CY 5 3 

CO M MO N /N03/N S , K SI ZB Y , N FLA , N -3S , MA XN , NF Pi 

C0MMCN/N06/NRBD, NBIFC, MXSMI, NFLAP, KPLAP, NOT, MXT2P1 ,NCSE 

COMMON/FKNl /FK'N 

COMMON/I PH1/NNFAZ, NNBS 

COMMON /EETS/H, FAH, FLH ,CTH 

COM MON/I SMB/ MS 

COMMON/AKTAU/AKCI (4) ,TA!J(4) ,SMLA (4 ) , A K (4 ) , AC (4) , CAPS, C? PC 
COM M0N/N05/NP,NSP, MO DS , M FP , ROUT ,MXQ , HAS , MBS , NET , N PS 
COM MON/ EN AM K/CS (4,2*0) ,SN (4, 30) 

EATA ALP/1,3,5,6,9,10/ 

JJ=NP+NFP1 
DO 60 1= 1 , MXS MI 
NMK=J J-I 
IL=I ABS { N M K ) 

IF(NMK) 16,15,17 

15 .EXP OF = (1 . ,0. ) 

GO TO 18 

16 A=CS (MS ,IL) 

E=SN (MS, IL) 

EXPOF=CMPLX (A, -B) 

GO TO 18 

17 A=CS (MS , NMK) 

B-SN (MS , NMK) 

EXPOF=C MPLX (A, B) 

18 CONTINUE 

KT1 = NBI FC* (1-1 ) 

IF (MS.GT.1) GO TO 35 
DO 30 JP = 1 ,-MXT2Pl 
30 FKN (KTl + JP) = (C . C , 0 . 0 ) 

35 KT = KTl + KXT2P1 + (MS-1) *NHBD 
DO 50 P = 1 ,6 
L = (1-1) *12 + ALP (P) 
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50 FKN (KT+P) =-H (L) *EXPOF 
FKN (KT+7) =-FAH (I) *EXPOF 

FKN (KT + 8 ) =SMLA (MS ) * ZKM (I , MS) *CTH(I) *EXPOF 
IF (MFLAP) 60,60,55 
55 FKN (KT + 9)=-FLH (I)*EXPOF 
60 CONTINUE 
BETURN 
END 
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CCMPLEX FUNCTION ZKM(I,MS) 

INTEGER CYl 8A 
INTEGER CY40,CY41,CY42 
INTEGER CY43,CY45,CY3,CY18,CY44 
COMPLEX CYl ,CY2,CY5,CY6 

COMMON/C¥I/CY3,CY18,C¥40,CY41,CY42,CY43,CY44,CY45,CY18A 
COMMON/CYC /CYl ,CY2,CY5,€Y6 

COMMON/C XR/CY4,C¥7,C¥8*CY9,CY10,CY11,CY12,CY13,CY14,CY15,CY16 # 

1 CY17,CY19,CY20,CY21,CY22,CY23,CY24,CY25,CY26,CY27,CY28,CY29, 

2 CY30,CY31,CY32,CY33,CY34,CY35,CY36,CY37,CY38,CY39 f CY46 # CY47,CY48, 

3 CY49#CY50,CY51,CY52,CY53 

COMMON /N03/NS ,NSIZBY ,NF1A» NES , MAXN, NFP 1 

COMMON/AKTAU/AKCI (4) ,TAU<4) , SMLA (4) / AK(4) ,AC(4) ,CAPK,CAPC 

KSM1=I-NFP7 

CY5=CYl-CY2*KSML*CY4 

ZKM=SMLA (MS) * (1+€Y5*TAU (MS) ) 

RETURN 

END 
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CVEELiY (NS 6, 1,6) 

PROGRAM Q5CLN 
COMPLEX DETSV 
COMPLEX FKN (12 3) 

COMPLEX QU {1 23) 

COMPLEX TEN (1681) 

COMMON/EPSA/QU, DETSV 
CCMMON/FKNl/IKN 

COMMON /N03/NS ,NSIZEY ,NFEA , NFS, MAXN,NFPl 
COM MCN/N04/N COLS , N B , NF , NEIFC,NBISC, MXTKN , NIG 
COM M0N/N05/NP,NSP, MODE , MFP , NOUT , MXQ , NAS, NBS , NET, NPS 
COM MCN/N06/NRB D, NR IFC , MXSMI, NFL AP, MFLAP , NCT , MXT 2P1 , NCSB 
COM M0N/NO7/I , J 
COM MO H/T KNl /TKN 
COMMON/DIMS/NN1 ,NN2,NN3 
C 

REWIND 1 

JJ = NP + NF Pi 

NCMAT = NR IFC* MXSMI 

WRITE (1) MXSMI, NRIFC, NF 

do io i = i , mxsmi 
DO 4 JP = 1 , NE IFC 
LL = (I-1)*HRIFC + JP 

4 QO (LL) = FKN (LL) 

IF (NN3 . EQ. 0) GO TO 3 
WRITE (6, 20OC) <00 (ID ,11=1, MXQ) 

3 CONTINUE 

DO 10 J = 1 , MXSMI 
CALL TKNS 

IF (NN3 . EQ. 0) GO TO 1C 

5 WRITE (6,2000) (TKN (II) ,11 = 1 , MXTKN) 

20O0 FORMAT (5X,8E14.6) 

10 WRITE (1) (TKN ( L ),L=1, MXTKN) 

1003 WRITE (1 ) (QU (L) ,L=1 , MXQ) 

1004 REWIND 1 

WRITE (6 , 2000) (QU (L) ,L=1,MXQ) 

2222 FORMAT (1H ,2E16,7) 

RETURN 

END 
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SUBROUTINES TKNS 
COMPLEX TKN (1681) . 

COMMON /IPH1/NNFAZ,NNBS 

COMMON /H04/ NCCLS , N B , N F , NEI PC, NEISC, MXTKN , NIG 

COMMON /N03/NS ,NSI ZE Y , N FEA , NES , MAIN , NFPl 

COMMON /NOS / N P , NS P, MODE, KFP , HOOT , MXQ , S AS , N3S , N ST ,N?S 

COM MCN/N06/NRBD, NRIFC , M XSMI , NFLAE, MFLAP, NCT, MXT2P1 , NC 

COMMON /NQ7 / I,J 

COMMON /N08/J YE 0 

COMMON /TKNl/ TKN 

DO 5 L=1 , MXTKN 
5 TKN (L) =(0. f 0.) 

IF (NSP) 10,10,25 
10 CALL BLA 
CALL BLB 
GO TO 60 

25 IF (JYEO. SQ. 1) GO TO 26 
CALL SNA 

GO TO 30 

26 CALL GYA 

30 IF (NNFAZ.FQ.O) GO TO 3 3 
IF (I.NE.KFP1) GO TO 32 

33 IF (I*NE.J) GO TO 32 
IF (JYTO.BQ. 1) GO TO 34 
CALL SWB 

GO TO 32 

34 CALL GYB 
32 CALL BLA 

CALL BLB 

IF (NB-1) 60,60 ,31 

31 CALL BLP 
60 CONTINUE 

BET URN 
END 



SUBROUTINE BLA 
INTEGER ALP (6) 

INTEGER Q 

COMPLEX B (1 8C0 ) ,SMLB (300) ,SMLC (300) ,SMLD (300) 

COMPLEX TK N (1681) 

COMMON /BTS 1/ B,S M LB , S MLC, SMLD 
COMMON /TKN1 / TKN 

COMMON /N02/ M XCSB ,NBSBC , MXCPM , NEBC 

COMMON /N03/NS , NSIZEY , NFEA , NES , MAXN , NFP 1 

COMMON /N04/ NCOLS , NB ,N F ,NBIFC , NEISC, MXTKN , NIG 

COM MCN/N06/NBBD, NRIFC, MXSMI, NFLAP, MFL1P , NCT, MXT.2P 1 , NCSB 

COMMON /NO 7/ I,J 

BATA ALP /I, 3,5,6,9,10/ 

DO 1§ KK=1,6 
DO 10 Q=1, NCOLS 

KT=N£IFC+HXT2P1+ <Q-1)*HRIFC + KK 

1= (J-1) *NEBC+ (1-1) *MXCPM+(Q“1)*12+ALP (KK) 

10 TKN (KT) =B (L) 

DO 15 KK=1 , 6 

KT= NEIFC +NC0LS*NRIFC+MXT2P1 +KK 
1= (J-1) *NESBC + <I-1) * MXCSB+ ALP (KK) 

15 TKN (KT) = -SMLC(L) 

DO 16 KK=1 ,6 

KT=NEIFC+ (NCOLS + 1 ) *NRIFC+MXT2Pl + KK 
L= (J-1) *NESBC+ (1-1) *MXC5B+ALP{KK) 

16 TKN (KT) =-SHLB(L) 

IF (MFLAP) 23/23,21 

21 DO 22 KK=1,6 

K1=NEIFC+ (NCOLS + 2) *NRIF C + MXT2P1 +KK 
L= (J-1) *NESBC + (I- 1 ) =* M X C5 B + A LP ( K K) 

22 TKN (KT) =-SMLD(L) 

23 CONTINUE 
RETURN 
END 
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SUBROUTINE BLB 
INTEGER Q 

COMPLEX FAB (150) ,FLB (150) ,CTB (1 5f>) 

COMPLEX FAS B (2 5) ,CTSB(25) ,FLSB{25) , FLSC(25) , 

1 FASD (25) ,CTSD (25) , FLSD (25) ,FASC (25) 

COMPLEX TKN (1681) 

COMPLEX YKM 

COMMON/AKTAU/AKCI (4) ,TAU(4) ,SMLA(4) ,AK{4) , AC{4) ,€APK,CAPC 
C0MM0N/BTS2/FAB # FL8,CTB 

COMMON/BTS3/F ASB ,C:TS B , FLS B , FLSC , FASD,CTSD, FASC 

COMMON /N03/NS,NSIZEY,NFEA,NBS, MAXN,NFP1 

CCMMCN/N04/NC0LS, NB , N F, NSIFC, NSISC , MXTKN,NIG 

COM MGN/NQ6/NRBD, NR IFC ,MXSMI,NFLAP, MFLAP , NCT , MXT2P1 , NCS B 

C0MM0N/N07/I , J 

COM MQN/TKN1/TKN 

MS = 1 

DO 10 Q=1, NCOLS 

KT= NSIFC +M XT 2P'5 + (Q-1) * NR IFC 

1= (J-1) *NCOLS*MXSMI+ (1-1) *N COLS +Q 

TKN (KT + 7 ) =FAB <L) 

TKN (KT+ 8) =- SML A (MS) *CTB (L) *YKM(I,MS) 

IF (MFLAP. EQ.O) GO TO 10 
TKN (KT + 9) =FLB (L) 

10 CONTINUE 

KT=NBIFC +NRIFC* (NCOLS + 1) +MXI2P1 + NCOLS 
L= (J-1) *MXSMI + (1-1) +1 
TKN (KT + 1 ) =-F AS B (L) 

TKN (KT + 2) =AKCI (MS) *CTSB (L) 

KTl =KT+NRIFC 
IF ( MFLAP . E0 • 0) GO TO 11 
IF (NFLAP.GT, NFEA) GO TO 12 
TKN (KTl + 1) =-FASD (L) 

12 IF (NFLAP.GT. NCT) GO TO 11 

TKN (KTl +2) = S ML A (MS) * YKM (I, MS) *CTSD (L) 

11 CONTINUE 

IF(MFLAP.EQ.C) GO TO 20 
IF (NFLAP.LE. NFEA) GO TO 14 
KT2 =NEI FC + M XT 2 P 1 + N HI FC* NCO L S+ NC CL S 
TKN (KT2 + 3) =-FLSC (L) 

14 IF (NFLAP.LK.NCT) GO TO 20 
TKN (KT+3) =- FLSB (L) 

20 CONTINUE 
RETURN 
END 
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SUBROUTINE BLP 
INTEGER P,Q 
COMPLEX EXPON 
COMPLEX TKN (16 81) 

COMPLEX YKM 

COMPLEX FAB (15 ) , FLB (150) ,CTB(15A) 

COMPLEX F AS B (25) ,CTSB (25) , FLSR (25) ,FLSC (25) , 

1 FASD (25) ,CT3D<25) ,FLSD{25) ,FASC(25) 

COMMON /AKTAU/ AKCI (4) ,TAU (4) ,3MLA (4) , AK (4) ,AC(4) ,CAPK,CAPC 
COM M0N/BTS2/FAB ,FL B , CTB 

CO M MO N/BTS 3 /F AS 3 , CTS B , FLS R , FLS C , F A S D , CT S D , F A SC 

COMMON /NO 3/NS ,N SI ZEY , NFBA , NES , MA XN,NFPl 

CCKMCN/N04/NCOL5,NB,NF,NEIFC,NiSISC,MX?KN,NIG 

CO M MO N /N06 /N P B D , N R IFC , M X S M I , NF L A P , MF L A P , NC T , MXT 2 P 1 , N CS B 

COMMON /N07 / I,J 

COMMON /TKN 1 / TKN 

NMK = J - I 

Kl = NEIFC + MXT2P1 

K2 = Kl + B 

DO 20 MS = 2, MB 

K3 = (MS-1) * (NBISC + ME 3D) 

DO 5 0 = 1 , N B B D 

K4 = Kl + (Q - 1 ) * NR I PC 

DO 5 P = 1, NRBD 

KT = K4 + P 

KT2 = KT + K3 

5 T KN (KT2 ) = TKN (KT) ♦ SXFON (NMK, MS) 

CO 10 Q = 1, NCOLS 
KT = K2 + K3 + (Q - 1) X NRIFC 
L*(J-1)*'BC0L5*HXSHI + (1-1) * NCCLS + Q 
10 TKN (KT) =-S ML A (MS) ♦CTB (L) ♦YKM (I, MS) ♦EXPON (NMK, MS) 

KT-K2 + K3+ (NCOLS + 1 ) ♦NRIFC 
L= (J-1) *MXS MI + (1-1) +1 

TKN (KT) = AKCI (MS)*CTSB(L) *BXPON (NMK, MS) 

IF(HFLAP.SQ.G) GO TO 20 
IF (NFLAP, GT. NOT) GO TO 20 
KTl =KT+ NRIFC 

TKN (KTl) = SKLA (MS) ♦YKM (I, MS) *CTSD (L) ♦EXPON (NMK, MS) 

20 CONTINUE 
RETURN 
END 
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SUBROUTINE SWA 
INTEGER P, Q , OS 
COMPLEX TKR (1681 ) 

COMPLEX XNLQ 
COMPLEX YKM 
COMPLEX ZLN 
COMPLEX VLN 
COMPLEX WNLQ 
COMPLEX EXPON 

COMMON/AKTA U/1KCI (4) ,TAU (4) f SMLA (4) ,AK(4) ,AC (4) ,CAPK,CAPC 

COMMON /N03/NS , NSI ZEY , NFE A , MBS, MAXN , NFP 1 

COMMON /R04/NC0LS , N B , H F, N'SIFC , HE ISC , MXTK N ,.NI G 

COM MQN/N06/NEBD # NB 1FC ,MX5MI , NFLAP # MFLAP # NCT , MXT2P1 ,NCSB 

COM M0N/NO7/1 , J 

COMMON/TKNl/TKN 

COM MON/I PHI /NN FAZ , NNBS 

COM MO N/AEPPM/DMS (4) 

COM MO N/SWASH/S WG J , SW El , S WM , S MB 
IF (I.NIS.J) GO TO 30 

THE I=J TEEMS ARE THE SMALL K=S MALL N TERMS 
DO 17 L=1 , M XT2P1 
LS = L-MAXN-1 
LL = (L-1) *NBIFC + L 
17 TKN(LL) = ZLN (IS, I) + VLN (I,LS) 

DO 20 MS = 1 , NS 
DO 20 L = 1 , MXT2P1 
LS = L-MAXN-1 

CFDL=1. + DMS (MS)* (1 , + (LS*LS- t) /(I . +LS*LS*S WGJ/SWEI) ) /SWR 
LL = (L-1) *N RIFC +MXT2P1 + (MS-1)*BRBD +8 
20 TKN(LL) = Y'K M (I, MS) *EXPON (LS , MS ) *CFDL 
GO TO 50 
30 IMJ = I-J 

DO- 18 L= 1 , MXT2P1 
DO 18 Q = 1 t MXT2P1 
IF (L.SQ.Q) GO TO 13 
LMQ = L-Q 

IF (IMJ, NE« LMQ) GO TO 18 
LS = L-MAXN-1 
QS = Q-MAXN-1 
LL = (L-1) *NRIFC + Q 

TKN(LL) = XNLQ (I,LS, QS) + W NLQ (I , I,S , QS) 

18 CONTINUE 
50 EE TURN 

BHD 



SUBROUTINE SWB 
INTEGER P,Q,QS 
COMPLEX TKN (1681) 

COMPLEX EXPON 

COMMON /N03/NS , NSI'ZEY ,NFEA , NES, MAXN, NFPl 
CO M MON /N04/NCQLS , N B , N F , NEIFC f N1 ISC , M XTKN , NIG 
CCMM0N/N05/NP, NSP, MODE, MFP, NOUT ,MXQ,NAS , NBS , NET , NPS 
COM M0N/N06/NRBD, NRIFC,MXSMI,NFLAP, MFLAP , NCT , MXT2P1 , NCS B 
COM MC N/N07/I , J 
COMMON/TKN1/TKN 

COM MON/AKT'AU/AKCI (4) ,TAU (4) ,SMLA<4) ,AK(4) ,AC (4) ,CAPK,CAPC 
COM MON/AERDM/DMS (4) 

COMMON/SWASH/S WGJ, SWEI , SWM, SWR 

IF (NBS.EQ.O) -GO TO 1 3 

KSML=I-NFP1 

NPMK=NP-KSKL 

DO 14 Q=1 , MXT2P1 

QS=Q-M AXN-1 

NMKQ=NPMK-QS 

IF (NMKO. EQ.S) GO TO 5 

RNMKQ=1.0*NMKQ 

RFA=RNMKQ/NBS 

NFA=N MKQ/NBS 

DIF =RF A-1 . u *NFA 

IF (DIF. GE. 0.0) GO TO 2 

DIF=-DIF 

2 IF (DIF. GT, .05) GO TO 14 
5 LL=N2IFC +NRIFC* (NCOLS + NCSB) + Q 

DFQR=1 . 0 + D.MS (1 ) * {1 + <QS*QS-1) /( ? +QS*QS*S WGJ/SWEI) ) /SWR 
TKN (LL) =-NBS*DFQR/SMLA (1) 

14 CONTINUE 
GO TO 23 

13 DO 21 M S = 1 , N B 
DO 21 Q=1 , NCSB 
DO 21 P=1 , MXT2P1 
QS= P-M AXN-1 

LL=NEI FC + (MS- 1 ) *UEISC +N3IFC* (NCOLS + NCSB) + (Q-1 ) *NRIFC+P 
DFQR = 1 + DMS (MS) *(1 + (QS*Q3~1) / (1 +QS*QS*SWGJ/5WEI) ) /SWR 
21 TKN (LL) = - BXP ON (-QS , MS) *DFQR/5 MLA ( MS) 

23 RETURN 
END 
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o o o o o n 


SUBROUTINE GYA 
INTEGER Q,QS 
COMPLEX YKM 
COMPLEX TKN <1681) 

COMPLEX GZLN 

COMPLEX GXN 

COMPLEX GXNK 

COMPLEX CY1,CY2,C¥5,CY6 

COMPLEX EXPON 

COM MON /NO 3/ N S , N SI ZE Y , N FEA , NES , M A XN , NFPl 

CO MMCN/N04/NC0LS , NB , N F, MEIFC, NEI5C, MXTKN , NIG 

C0MM0N/N06/NRBB, NF.IFC , MXSMI, NFLAP, MFLAP , NCT , MXT2P1 , NCSB 

COM M0N/N07 /I , J 

COMMON/TKN1/TKN 

COMMGN/IPHi/NNFAZ, NNBS 

COM MON/GYR/GYM ,GYK,GYC,GKB, GCB , GKP , GCP , GRP,GIX,GIY,GIZ 
COM MON/C YC/C il ,CY2,CY5,CY6 

COMMON/C YR/CY4 , CY7 ,C Y8 ,CY9 ,CY10 ,CYl 1 , CY12 ,C Yl 3 , CY14 , CYl 5, CYl 6 , 

1 CYl7,CYl9,CY20,CY21,C¥22,CY23, CY24, CY25,CY25, C Y27,CY28 ,C Y29 , 

2 CY30,Cy3l,CY32,CY33,CY34,CY35,CY36,CY37,CY38,CY39,CY46,CY47,CY48, 

3 CY49,CY50,CY51,CY52,CY53 
COMMON /RNAME/ CS (4 ,20) ,SN (4 , 20) 

JMI=J-I 

IF (JMI.NE.O) GO TO 30 


COMPUTES DIAGONAL ELEMENTS OF DIAGONAL BLOCKS 

KSML=I-NFP1 

C Y5=CY1“CY2*KS ML i *CY4 

CY6=CY5*CY5 

DO 5 L=1,MXT2P1 

LS=L-MAXN-1 

LL= (L-1) *NRIFC+L 

IF (LS) 2,3,4 

2 GZLN=CY6*GIX+CY5*GCB+GK3+CY7* (GIZ-GIY) 

GO TO 5 

3 GZLN=CY6*GYM+CY5*GYC+GYK 
GO TO 5 

4 GZLN=CY6*GIY+CY5*GCB+GKB+CY7* (GIZ-GIX) 

5 TKN (LL) =GZLN 

COMPUTES OFF DIAGONAL ELEMENTS OF DIAGONAL BLOCKS 

IF (MAXN.EQ. 0) GO TO 10 
Q=MXT2P1 
L = M X T 2 P 1 

TKN (Q) =-CY4* (CY5* (GIZ-GIX-GIY) -GCB) 

LL= CL- 1) *NRIFC+1 
TKN (LL) = -TKN (Q) 
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o o O on 


CO MPUTLS BLADE 


roaPLING TEEMS 


10 CONTINUE 

DO 2D KS=1 # HB 
DO 20 L=1 # SXT2P1 
LS=L-'HAXN-1 

LL= (L-1) *NEIFC+HXT2P1 + (MS-1 ) *HRBD+8 
IF(LS) 6,7,8 

6 TKN (LL) =-GRP*SN (MS,1)*YKM (I, MS) 

GO TO 20 

7 TKN (LL) =-YKM (I,MS) 

GO TO 20 

8 TO (LL) =GPF*CS (MS , 1 ) * YK M (I , M S) 

20 CONTINUE 

GO TO 50 


COMPUTES OFF DIAGONAL TERMS IN OFF DIAGONAL BLOCKS 

30 IF ( (JHI-2) .£.0.0) GO TO 35 
IF { (JMI+2) • 3Q. 0) GO TO 35 
GO TO 50 
3 5 M X= MA XN+ 1 
KSML=J-NFPl 
CY5 = CY1-CY2*KSML*CY4 
GXNK-GKP+CY5*GCP+.5*CY2*CY4*GCP*JMI 
DO 40 L=1 , MXT2 Pi 
DO 40 Q=1, MX T 2 Pi 
IF (L. EQ. MX) GO TO 40 
IF (Q.FQ.MX) GO TO 40 
LS-L- MAXN-1 
Q£=Q-MAXM- 1 
I L= (L-1 ) *NRIFC+Q 

40 TKN (IL) =-. 25*GXNK* (LS+QS-C Y2*LS* (LS-QS) *JMI*.5) 

50 CONTINUE 
RETURN 
END 
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SUBROUTINE GYB 
C 

C CANNOT HANDLE N BS GE EATER THAN ZERO 

C 

IHTEGFB Q,QS 
COMPLEX TKN (1681) 

COM MGN/N03/N S, NSIZEY ,NFEA , NHS, MA XN,NfPl 
C0 HM.cn/B 04/HCOLS # NB # N E , NBIFC, NEI3C, MXTKN , NIG 
COM HON/NOS /BRED, KB IFC,MXSMI # NFLAP, MFLAP , NCT , MXT2P1 , NCS B 
CCMMCN/N07/I , J 
COM MON/IK Ml /TKN 

COHMOA’/AKTSU/AKCI (4) ,TAU(4) , SMLA(4) ,AK(4) ,AC(4) ,CAPK,CAPC 
COMMON/BNAMB/CS (4,20) ,£N (4,23) 

CO M HON/G YR/GYK ,GYK,GYC,GKB, GCB , GKP , GCP, GRP , GIX, GI Y, GIZ 
C 

DO 20 MS = 1, NB 
CO 20 0=1 ,MXT2Pl 
QS=Q- MAXN-1 

LL=HEIPC + (MS-1) *NEISC + NEIFC« t (NCOLS+1) + Q 
IF (QS) 2,4,6 

2 TKN (LI) =GF P* SN (MS,1) /SHLA (MS) 

GO TO 20 

4 TKN (LL) =1.D/SMLA (MS) 

GO TO 20 

6 TKN (I L) =-GSP*CS (MS ,1) /SHLA (MS) 

20 CONTI NOE 
RETURN 
END 
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COMPLEX FUNCTION EXCHI (L,Q,J) 

C CHE ATE EXP (I* (L-Q) *CHI { J) ) 

INTEGER 0 

COM MON/E NAM El/CSl (4,06) ,SNl(4,06) 

LQ=L-Q 

ILQ=I A3S (LQ) 

IF(LQ) 16,15,17 

15 BXCHI=(1. ,0.) 

GO TO 18 

16 A= CS1 (J,ILQ) 

B-SNl ( J , XL Q ) 

BXCHI = CMPLX (A,-B) 

GO TO 18 

17 A=CS1 ( J , ILQ) 

B=SN1(J,IL0) 

EXCHI =CMPLX (A, B) 

18 CONTINUE 
RETURN 
END 
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COMPLEX FUNCTION EXPON(L,MS) 

CREATE EXP (I*L*PHIM) 

COMMON/RNAME/CS (4,20) , SN (4,20) 
II=IABS (L) 

IF (L) 16,15,17 

15 EXPON= (1, , 0 , ) 

GO TO 18 

16 A=CS (MS , IL) 

B = SN (MS, IL) 

EXPON=CMPLX (A, ~B) 

GO TO 18 

17 A=CS (MS,L) 

B = SN ( MS , L) 

EXPON=CMPLX (A, 8) 

18 CONTINUE 
RETURN 
INI) 



COMPLEX FUNCTION VLN{I,LS) 

COMPLEX VN , V Nl , VN2 

INTEGER CYl 3 A 

INTEGER CY 40 , C Y4 1 , C Y 4 2 

INTEGER CY43,CY45,CY3,CY13, CY44 

COMPLEX CYl ,CY2,CY5,CY6 

COMMON /AERBP/ BJ (4) 

COMMON /ABB TP/ AKT (4) f ACT (4) ,AKP (4) ,&CP (4) 

COM MON/AKTAU/ AKCI (4) , T A r J (4) , SML A {4} ,AK (4) ,AC (4) ,CAPK,CA PC 
COM MON/C YI/C Y3 ,€ Y‘i 8 , CY40 , CY 4 1 , CY42, CY43 , CY4 4 , CY 45 ,CYl 8A 
COM HO N /C YC /C Y 1 , C Y2 , C Y 5 , C Y 6 

COM HON/C YB/CY4,CY7,CY8,CY9,CYltf,CY?1,CYl2 f CYl3,CYl4,CYl5,CYl6, 

. 1 CYl7,CYl9,CY2v,CY;>1,CY22,CY23,CY?.4 # CY25,CY26,CY27,CY28,CY29, 

2 CY30 f CY31 ,CY32.,CY33,CY3 4, CY35, CY36,CY37 ,CY38,CY39,CY46 ,CY47,CY48, 

3 CY49,CY50,C¥51,CY52,CY53 

COMMON /N03/NS ,NSIZBY,NFEA, IBS , MAIN , NFPi 

COMMON /SWASH/ SWG J, SWX1 , SWM ,SWR 

B=SWR 

KSML=I-NFPl 
C Y5 = C Y1-CY2 * KS ML*CY4 

CFL = 1 . + (LS*LS“ 1 ) / (1 + LS*LS*S WGJ/S WE I ) 

VN= (C . ,C.) 

VN1 = (0. ,0.) 

V'N2= ( 0 . ,0.) 

UO 5 J=1 ,NES 

VN=VN + AKT (J) +(CY5-CY2*LS*CY4) ♦ACT (J) 

V N 1 = V N 1 + ( AK ( J) + (CY5-CY2*LS*CY4) *AC (J) )*B.J («T ) * (R-BJ (J) *CFL) 

5 VN2-VN2+AKP (J) + (CY5-CY2*LS*CY4) *AC? (J) 

VLB = (L$*1$*VH-VN1*CFL+CFL*CFL*VN2)./(R*R) 

BET URN 
END 
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COMPLEX FUNCTION W NLO {I , IS , QS) 

COMPLEX ‘fiXC HI 

COMPLEX WS ,WN1 ,WH2 

INTEGER QS 

INTEGER CY1 8 A 

INTEGER CY40,CY41,CY42 

INTEGER CYt*3,CY45, CY3 , CYl 8 , CY44 

COMPLEX CY1 ,CY2,CY5,CY6 

COMMON /ABB BP/ BJ(4) 

COMMON /ABB TP/ AKT {4 ) , ACT (4) , A K P (4) , ACP (4 ) 

COMMON/AKTAU/AKCI (4) ,TAU(4) ,SMLA(4) ,AK{4) ,AC(4) ,CAPK,CAPC 
COMMON/CYI/CY3 ,C Yl 8 , C Y40 ,C Y 41 , C Y42 , C Y43 , CY4 4 , CY 4 5 , C Yl 3A 
COM MO N/C YC/CYl ,CY2,CY5,CY6 

COM MO N/C YP/CY4 , C Y7 , C Y3 , CY9 , C Y 1 1 , CY1 1 , CY 1 2 , C Yl 3 , CY 1 4 , C Yl 5 ,C Y 
t CYl7 f CYl9,CY20 # CY2'l ,CY22 ,C Y23 , CY24 , CY25 ,CY26 ,CY27 , CY2S ,CY2 

2 CY30,CY31 ,CY32,CY33,CY34, CY35, CY36 , CY37 ,CY38, C Y 3 9, CY46 ,CY4 

3 CY49,CY50,CY51,CY52,CY53 

C CM MO N /N03/NS , NSIZBY , NFEA , NBS , MAXN , N FP 1 
COMMON /S'JASH/ SWGJ ,SHEI ,SWM ,SWR 
F=S WR 

KSML=I- NFP1 

CY5 = CY1-CY2*KSML*CY4 

CFQ=1 . + (QS*QS-1) / (I . +QS*QS*SWGJ/SWi?I) 

CFL = 1. + .{LS*IS-1)/(1.+LS*I,S*SW6J/Sil3I) 

«N= (IK ,0.) 

w»1 = (0. ,0.) 

WN2=(C. ,0.) 

BO 5 J=1,NBS 

WN=WN+ (AKT (0) + (CY5-CY2*QS*CY4) * ACT (,J) ) *SXCH.I (LS ,QS, J) 

WN1 =WN1 + (AK (J) + .CCY5-CY2*QS*CY4) t-AC (J) )*BG (J) * (B-B J { ,J) *CFQ) 

1 *EXCEI (LS , QS , J) 

5 WN2 = WN2+ (AKP (J) + (CY5-CY2*QS*CY4) *ACP (J) ) *BXCHI (LS ,QS, J) 

WNLQ = <QS*LS*NN-WN1*CFL+CFQ*CFL*WN2)/{B*B) 

RETURN 

END 
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COMPLEX FUNCTION XNLQ (I , LS , QS) 

INTEGER QS 

COMPLEX EXCHI 

INTEGER CY18A 

INTEGER CY4G,CY41,CY42 

INTEGER CY43,CY45,CY3,CY18, CY44 

COMPLEX CY1 ,CY2,CY5,CY6,XN,XNl ,XN2,XN3 

COMPLEX C5 , 

COM MON/A1R BP/B J (4) 

COMMON/AKTAU/AKCI (4) ,TAU (4) , SMLA (4) , AK(4) ,AC (4) ,CAPK,CAPC 
CCMHON/CYI/C Y3 ,CYl 8,CY49,CY41 ,CY42,CY43 ,CY44 ,CY45 ,CY1 8A 
COMMON/C YC/CYl ,CY2 # CY5 # CY6 

COM MON/C YR/CY4,CY7, CY8 , CY9 , CYtQ , CY1 1 ,CY 1 2 ,CYl 3, CYl4,CYl 5, CYl 6 , 

1 CY17,CY19,CY2Q,CY21,CY22,CY23,CY24 ,CY25,CY26 ,CY27,CY28 # CY29, 

2 CY30,CY31 ,CY3-2,CY33,CY34,CY35, CY36 , CY3 7 , CY38,C Y39, CY46 ,CY47 # CY48, 

3 CY49,CY50,CY51,CY52,CY53 

COMMON /N03/NS,NSIZEY,NF£A,MES,MAXN,NFPl 

COMMGN/N04/N COLS, NB , N F, NEIFC , N3ISC, MXTKN, NIG 

COMMON /SWASH/ S WG J , S WSI , SWM , S W R 

KSML=I-NFP1 

C Y5 = C Yl-CY2*KS ML^CY4 

I F (QS . SQ, LS) GO TO 15 

XHLQ=(0./0.) 

CFQS - (1 . + (QS«QS-1) / {I. +QS*QS*S W GJ/S WEI) ) /SWR 
CFLR= (1 . +(LS^LS-1) /(I . +LS*L3*SWGJ/SWEI) ) /SWR 
DO 10 j=1 , NES 

10 XNLQ=XNLQ+ (AK (J) + (CY5-C Y2#Q3*CY4) *AC (J) ) *HXCHI (LS ,QS ,J) * 

1 (1 ,-BJ ( J) ♦CFOS) 

IF (NIG .SQ. 0) GO TO 16 
C 5 = C A P K + C A PC* ( C Y 5- C Y 2 * Q S '*C Y 4 ) 

XN1 = (0.,,0.) 

X N 2 = ( I) . f 0 . ) 

XN3 = (0. ) 

DO 12 J = 1 , N E S 

XN= AK (J) + (CY5-CY2*QS*CY4) *AC<J) 

XN 1 =XNl +XN 

X N2 = X N2 + X N* E XC HI ( LS , 0 , J) ♦ { 1 . - B J ( J) * CFLR ) 

12 XN3=XN3 + XN*EXCHI (0,QS, J) (J) *CFQR) 

X NL Q= X NL Q-XN3*XN2/(C5 + XN1) 

GO TO 16 

15 XNLQ =(■"). ,0 . ) 

16 CONTINUE 
RETURN 
END 
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COMPLEX FUNCTION YKM(I,MS) 

INTEGER CY18A 
INTEGER CY40,CY41,CY42 
INTEGER CY43,CY45,CY3,CYl8,CY44 
COMPLEX CY1 ,CY2,CY5,CY6 

COMMON/C YI/C Y3 ,C Y1 8,CY40,CY41 ,C Y42 , CY43 , C Y4 4 ,CY45,CYl8A 
CCMMON/CYC/CYl ,CY2,CY5 ,CYo 

COM MON/C YR/CY4,CY7,CY8,CY9,CY10,CY11 ,CYl 2 ,C Y1 3 , CY14 , CYl 5, CYl 6, 

1 CY1 7 r CY1 9 , CY20, CY21 ,CY2 2, CY23, CY2.4,CY2 5,CY26,C Y27,CY28 ,CY29 , 

2 CY30,CY31,CY32,eY33,CY34,CY35,CY36,€Y37,CY38,CY39,CY46,CY47 f CY4S, 

3 CY49,€Y50,CY51,CY52,CY53 

COMMON /N03/NS ,NSI ZEY ,NFEA , NSS # MAXN , NFPl 

COM MON/ AKTAU/ A KCI (4) ,TAU'(4) ,SMLA (4) , AK(4) , AC (4) ,C.APK,CAPC 

KSML=I-NFP1 

CY5=CY1~CY2*KSML*CY4 

Y KM=S ML A (MS) * (1+CY5*TAU (MS) ) 

RETURN 

END 
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COMPLEX FUNCTION ZLN(LS,I) 

INTEGER CYl 8 A 

INTEGER CY40,CY41 ,CY42 

INTEGER CY43,CY45,CY3,CY18,CY44 

COMPLEX CY1 ,CY2,CY5,CY6 

COMPLEX 2X CHI, CLUJ 

COMPLEX Cl ,C2,C5,C6,C7,C3,C9,C10 

COM MON/ AER BP/3 J (4) 

COMMOW/AKTA U/AKCI (4) ,TAU (4) ,SMLA(4) ,AK (4) ,AC (4) ,CS*>K,CAFC 
C 0 M M 0 N / C Y I / C Y 3 , C Y 1 8 , C Y 4 0 , C Y 4 1 , C Y 4 2 , C Y 4 3 , O' Y 4 4 , C Y 4 5 , C Y 1 8 A 
CCMMQN/C YC/CYl ,CY2,CY5,CY6 

COM MO N /C YE /C Y 4 , € Y 7 , C Y 8 , C Y 9 , C Y 1 0“ , C Y 1 1 , CYl 2 , C Yl 3 , CY14 , CY1 5, CYl 

1 CYl 7 ,CYi S , C Y2*> , GY21 , CY2 2, CY23, CY24, CY2 5, CY26 , C Y27 , C Y28 ,C Y29 

2 CY39 ,CY 31 ,CY3 2,CY33,CY34 ,CY35 ,CY36 ,CY37 ,CY3 8 ,CY3 9, CY46 ,CY47 

3 CY49,CY5%CY51 ,CY52,CY53 

COMMON /MO 3 /N S , N S 1 7,¥« Y , N FEE , NE S , MAIN, NFPl 
C0MMGN/N04/NC0LS , N3, N F, NEIF C, N EISC, MXTKN , NIG 
COMMON /SWASH/ S WG J , S WEI , SW M , SWB 


S = SWR 

KSML=I-NFP1 
CY5=CYl-CY2*KSML*CY4 
C Y6 =CY5*CY5 

Cl = SWM* <CY6-CY8*LS*CY2*CY5-LS*LS*CY7) 

C2= (0. ,0.) 

C9= ( \ ,0, ) 

CFLB= (1 . + (LS*LS-1) /(I + L S* : LS* S WG J/S WEI ) ) /SWH 
DO 10 J-1 , NES 

CLN J = AK (J) +CY5*AC (J) -CY2*LS*CY4*AC (J) 

C2 = C2 + ^LN J 

10 C9=C9 +CLNJ* 4 (1-BJ (J) *CFLP.) 

C10=C2 

C2-C9 

IF (NIG .BO. C ) GO TO 11 
C2-C10 

C5= C A PK+ C A PC* (CY5-CY2*L3*CY4) 

C6= (0. ,0.) 

C7= (0, ,0.) 

DO 8 J=1 , NFS 

C8=AK (J) +CY 5* AC (J) -C Y?.*LS*C Y4*AC (J) 
C6=C6*C8*FXCIII (LS,3, J) * (1. -BO (J) *CFLB) 

8 C7=C7+C8*EXCHI (0,LS,J)* (1-BJ{J) *CFLR) 
C2=C9-C7*C6/ (C5+C2) 

11 IF { M A X N , EQ. 1) GO TO 12 
KX=1-LS*LS 

C3=2. *CY12*LS*LS*KX*KX 
C4=B*K*K* n ,/SWGJ+LS*LS/SWEI) 

C3=C3/C4 

ZLN=C1+C2+CMPLX(C3,0. ) 

GO TO 13 



12 ZLN=C1+C2 

13 RETURN 
END 
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OVERLAY (NS6,2,0) 

PROGRAM SOLVE 

COMPLEX DB (123,123) , DTPHAS , DPIVOT 
CCMELEX B (123,123) 

COMPLEX EPS (123) 

COMPLEX DETSV 
EQUIVALENCE (DB,B) 

COMMON /EPSA/ EPS, DETSV 
COMMON/I W A 4/NN 4, NN 5, NN6 

COMMON /CDETRM/ DPIVOT , DTPHAS , DTLGl 0, IDET 
REWIND 1 

READ (1) MXSMI,NRIFC,NF 

NORDER=HXSHI*NRIFC 

CALL INRAY (B, MXSMI, NRIFC) 

READ (1) (EPS (I) ,1 = 1 ,N05DER) 

REWIND 1 
IDET=0 

N=N ORDER 

CALL DC MAT (DB,N,EPS) 

WRITE (6 , 700) (EPS (I) ,1 = 1 ,NORDER) , DETSV 
5 CONTINUE 

700 FORMAT (5X,8E14. 6) 

RETURN 

END 
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SUBROUTINE INRAY ( B, MXS MI, N RIFC) 

COMPLEX B (123,123) 

DO 1 L= 1 , MXS MI 
IRO«?=L*NEIFC 
ISOWS=IRO WF-NRIFC+1 
DO 1 K=1 , MXS MI 
ICO L F =K * NR I FC 
ICOLS=ICOLF-NBIFC+1 

1 READ (1 ) ( (B (I, J) ,I=IEOWS,IROWF) ,J= ICOLS,ICOLF) 

RETURN 
END 



SUBROUTINE DC MAT {A ,N , Y) 

COMPLEX A (123/123) ,Y (123) , DOST , CAIJ , AMX , DONS, DYI , TEMP , DAK J, 
1 DY K , E h K K , DA IK , DPI V OT , DPH AS 
COMPLEX X (123) ,!K,GNB 
DIMENSION ICHG (123) 

COMMON /CDETRM/ DPIVOT , DPHA3 , ADET, IDET 

NDI M=1 2 3 
ADET=C. 0 
DSI GN=1 . 0 

DPHAS=CMPLX (1.0,0. 0) 

G NE = 1 .0 

ZEE 0 = 0, C 

API = N 

IF(ID??.J5Q.O) GO TO 651 
HP1-N+1 

DO 651 1=1, NPl 
651 X (I) =A (NPl ,1) 

650 CONTINUE 

DO 118 K=1 , N 
AMX = A (K, K) 

I MX = K 

DO 100 I=K , N 

IF (CABS (A (I, K) ) . LE. CABS (AMX) ) GO TO 100 
AMX = A(I,K) 

I MX =1 

100 CONTINUE 

102 IF (I MX. EQ. K) GO TO 136 
DO 104 J = 1 , N P 1 
TEMP=A (K , J) 

A (K , J) =A (IMX , J) 

104 A (IMX, J) =TEMP 
ICHG (K) =IMX 
TEM F=Y (K) 

Y (K) = Y (IMX) 

Y (I MX) = TEMP 
DPHAS=-DPHAS 
GO TO 108 

1C 6 ICHG ( K) =K 
1C 8 CONTINUE 

D AKK=A (K , K) 

901 FOP SAT (1X,I5,2E40. 16/) 

AMAG=CABS (DARK) 

IF (AMAG.NE, ZERO) GO TO 6 
WRITE (6,7) 

7 FORMAT (*5 MATRIX IN DC MAT IS SINGULAR*) 

STOP 

6 CONTINUE 

ADST=ADET+ALOG10 (AMAG) 

D PH AS = D PU AS * E A K K/ A MAG 
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DONE=CMPLX (1 «C' ,0.3) 

eakk=done/dakk 
DO 110 J=1,NPl 
110 A<K,J) "A (K, J)*1>AKK 
ft (K , K) = D A K K 

IF (IDET. SO. 0 ) GO TO 652 
TK = X (K) 

DO 653 J=K,NPl 
653 X (J) =X(J) rT K*A <K,J) 

652 CONTINUE 
DYK=Y (K) 

Y (K) =DYK*DAKK 
DO 114 1=1, N 

IF (I.EO.K) GO TO 114 
DftIK=A (I, K) 

DO 112 J=1 , NPl 

112 A (I, J) =A (I, J) -I>AIK*A {K, J) 

CALL ROW S U K { NP 1 „ N D I M , A ( 1 , 1 ) , A (K ,1) , DAXK) 
A (I, K) =DAIK 
D YI = Y (I) 

DYK = Y (K) 

Y <I) =DYI-DAI K*DYK 
114 CONTINUE 

DO 116 I -1 , N 
116 A (I,K) =~A (I,K) *DAKK 
A (K,K) =DAKK 
118 CONTINUE 

DO 12 2 K=1 , N 
L=N + 1-K 
KI=ICHG (L) 

IF (L.BQ.KI) GO TO 122 
DO 120 1=1 , N 
TEMP = A .(1,1) 

A (I ,'L) = ft (I , KI) 

120 A (I, K I) = TEMP 
122 CONTINUE 

IF (ID 3T.NB. 0) DPIV0T = X (NPl ) 

124 RETURN 
END 
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Machine Compatibility 

The Swashplate Dynamic Response Program has been run on 
NASA-Langley ' s CDC 6600 and CDC 6400 systems. The program was 
developed utilizing standard FORTRAN IV and is also WATFIV 
compatible. The program has not been run in its present form on 
an IBM 360/65. To execute this program on an IBM system the 
overlay cards would have to be removed and replaced with an IBM 
overlay structure and the PROGRAM cards changed to SUBROUTINE 
cards (for example, the PROGRAM SETUP card would have to be re- 
placed with a SUBROUTINE SETUP card) . In addition, the core 
requirements may be restrictive, especially if it is necessary 
to convert the program to double precision for accuracy of result 
to be equivalent to that obtained on a CDC system. 
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